
(NAS ft -CP-157128) HEL1CGYR0 PEELIMINAFY 
DESIGN, PHASE 2 Final Report 
(MacNeal-Schvend ler Corp.) 192 p HC A09/NT 
ftOI CSCl 22E 


Unclas 

G3/15 20705 








HELIOGYRO PRELIMINARY DESIGN 


PHASE II FINAL REPORT 


19 January 1978 


MS404-5 


L-L-LA 

THE MACNEAL-SCHWENDLER CORPORATION 






HELIOGYRO PRELIMINARY DESIGN 
PHASE II FINAL REPORT 

MS404-S 19 January 1978 


MS404-5 


HE! IOGYRO PRELIMINARY DESIGN 
PHASE II FINAL REPORT 


for 


Jet Propulsion Laboratory 
Pasadena, California 


Under Contract 954709 


19 January 1978 


THE MACNEAL-SCHWENDLER CORPORATION 
7442 North Figueroa Street 
Los Angeles, California 90041 


HELIOGYRO PRELIMINARY DESIGN 
PHASE II FINAL REPORT 

INTRODUCTION 

This report summarizes the work performed under Mod. 5 of Contract 
954709 for the Jet Propulsion Laboratory. The work included several 
aspects of the preliminary design of a Heliogyro Solar Sail Module for 
the Halley Rendezvous Mission. Mr. William Ruff was the Technical Pro- 
ject Manager at JPL. The MacNeal-Schwendler Corporation (MSC) was the 
prime contractor and their work was directed by MSC's president, 

Dr. Richard H. MacNeal. Astro Research Corporation (Astro) was a 
principle subcontractor and their work was directed by Astro's president. 
Dr. John M. Hedgepeth. 

The project, which began in February 1977, wss based on the results 
of a short conceptual design study (contract 954680) conducted during 
December 1976 and January 1977 by MSC and Astro. The objective was to 
develop a Baseline Design of a Heliogyro Solar Sail Module on a short 
schedule, to be considered for the Halley Rendezvous Mission in competi- 
tion with a Square Solar Sail Module and, later, in competition with a 
Solar Electric Propulsion System. The main features of the Heliogyro 
Baseline Design, including the work done at JPL as well as that done at 
MSC and Astro, are reported in the "Solar Sail Technology Readiness 
Report," JPL Report 720-1, 18 July 1977. The main part of the work 
done by MSC and Astro is reported in MSC's report number MS404-1, "Helio- 
gyro Preliminary Design, Final Report," August 24, 1977. 

The present Phase II Final Report covers work done since that date 


on the following six topics: 


» 



1. Design and analysis of a stowable circular lattice batten for 
the Heliogyro blade. 

2. Design and analysis of a biaxially tensioned blade panel. 

3. Definition of a research program for micrometeoroid damage to 
tendons . 

4. A conceptual design for a flight test model of the Heliogyro. 

5. Definition of modifications to the NASTRAN computer program 
required to provide improved analysis of the Heliogyro. 

6. A User's Manual covering applications of NASTRAN to the Heliogyro. 

Separate memoranda and technical notes on these topics constitute the 
body of the Phase II Final Report. 
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INTRODUCTION 


The Heliogyro concept was created in the mid 1960 ' s by Astro 
Research Corporation (Astro) and MacNeal-Schwendler Corporation 
(MSC) under the support of NASA Headquarters. References 1 and 2 
summarize most of the results of that work. 

Astro has a subcontract with MSC to support JPL in applying 
the Heliogyro concept to the Halley Comet Solar Sailing Mission. 
William Ruff is the Technical Project Manager at JPL. The work at 
MSC is being led by Richard H. MacNeal - Program Manager, and at 
Astro under the direction of Karl Knapp - Program Manager. John M. 
Hedgepeth is Astro's Senior Scientist. 

There are 12 blades in the Heliogyro design, and each blade is 
envisioned to be 8 meters in width and 7,500 meters in length. The 
blades are expected to be composed primarily of a thin membrane 
constructed of material such as Kapton film with an aluminum reflec- 
tive coating on one side and an infrared emissive coating on the 
other. An overall view of the baseline Heliogyro is shown in 
Figure 1. 

This report is one of a series dealing with the design and 
fabrication concepts of Heliogyro blades. In particular, this 
report details the Heliogyro stowable circular lattice batten. 
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DESIGN 


The Heliogyro blade requires chordwise battens at about 75- 
meter intervals to hold the leading- and trailing-edge tension 
members apart. These battens must be 

1. capable of resisting a design limit compressive load of 
approximately 6 newtons and temperatures of 600 kelvin, 

2. tightly stowable on the blade reel, and 

3. self-deploying when the blade is unrolled. 

In addition, the batten should be of lattice construction so that 
the thermal gradients through the cross section are small and so 
that the reradiation from the batten will not produce hot spots in 
the blade reflecting membrane. 

The batten design that meets these requirements is shown in 
Figure 2. I' is constructed of thin graphite/poly imide rods in a 
circular, cylindrical lattice. There are six square longerons that 
support the compressive load, and six 30-degree helical half-round 
spirals, three in each di~ection, that provide support to the lon- 
gerons to avoid local buckling. The longeron-spiral and spiral- 
spiral intersections are staggered to aid fabrication (by avoiding 
the criss-cross pileup) and minimize the local unsupported length 
of the longeron. 

The batten is fabricated in two halves and hinged together on 
assembly. The hinges allow full compaction when stowed and expan- 
sion to a circle when unrolled from the stowage reel. 

Two important principles are included in the design: 

1. Each half-batten is fabricated on a mandrel with a smaller 
diameter than that of the finished batten. 





2. The "hinges" are composed of interlacing fingers with no 

hinge pin required. 

These two principles allow the batten, when fully deployed, to 
behave structurally as if there were no hinge and the spirals were 
continuous. Furthermore, they permit fabrication without the neces- 
sity of complicated joints. The tendency of each half-batten to 
reach a diameter smaller than the batten diameter causes a preload 
on the hinges, permitting moment carry through. It also produces 
a force which actually holds the two halves together, resisting 
possible spreading forces. 

These principles are not new. They have been discovered and 
developed previously by Astro Research Corporation for solid tubular 
booms such as the Tablock BI-STF.M. The application of the principles 
to a lattice configuration is a new and important development. 

Figure 3 shows a sample batten segment fabricated as a part of 
the development of the design. It is constructed of graphite/ epoxy 
(graphite/poly imide roving was not available) arid demonstrates the 
geometry of the lattice halves and finger hinges. 

An accurately constructed model which was used to demonstrate 
the strength of the batten design for local buckling is shown in 
Figure A . This model was fabricated of graphite/epoxy as a complete 
circle with no hinges. It is 12 cm in diameter and consists of: 

longerons 0.48-mm square 

spirals 0.96-cm diameter half-round 

The calculated buckling strength of the longerons cu'tsxcered 
to be simply-supported columns between the intersections is 

P 4 

theor 


24.8 N 


An estimate of the general- instability buckling strength 
(involving deflection of the spirals) can be obtained by starting 
with Reference 3 and substituting in the appropriate stiffness of 
the longerons and spirals. in this process, the axial stiffness of 
the spirals must be ignored since there is no skin to cause cylinder 
action (mid-plane stretching) to occur. The theoretical general- 
instability load thus obtained after some effort is 

P = 20.6 N 

GI 

The measured strength is 

P = 22.6 N 

meas 

which is very close to the theoretical strengths. The observed 
buckling mode involved cross-sectional deformation. The results 
of the test establish an experimental determination of the amount 
of spiral stiffness required to make the general-instability load 
approximately equal to the local-buckling strength. 
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ANALYSIS 


Let longerons be square with thickness t^ . Let the helix 
angle be 30 degrees and the ratio of stiffener area to longeron 
area be r . The buckling loads are then 


Euler column 
buckling: 


,2 2 2 

3tt ER t^ 


EU 


(1) 


Local longeron 
buckling: 



( 2 ) 


Equation (2) is derived by assuming that the longerons go through 
the spiral intersections. By staggering the longerons and inter- 
sections, we get a potential local buckling load four times this 
much. Letting P./P = 1 actually gives a factor of four on local 
buckling provided that enough stiffness is incorporated in the 
spirals to produce stabilization of the intersections and avoid 
general instability. 

Manipulating Eqs. (1) and (2) gives 




2 2 

81 n e 6 

8 2 
L 


(3) 



2 2 I 

^ 81tt E j 


(4) 
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3tt e t 

“V 



The volume of material in the longerons and spirals is 


V 


6L t, 


(1+2r A» 


V = 4(1+ 2r A ) (3rr 2 E 2 )" 1/3 L 5/3 (P t P E]J ) 1/3 

Let 

L = 8 m 

E = 124.8 x 10 9 N/m 2 

P = 1550 kg./m 3 

V = 1.656 x 10 _6 (1+ 2r A ) (P^P £u ) 1/3 m 3 

t^ = 1.858 x 10 _4 (P^P eu ) 1/6 m 

R = 2.24 x 10" 2 (P^)" 1/6 m 
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Mass = 0.00257 (1 + 2r )(P,P ) 1//3 kg 

A t EU 


(ID 


where P^ and P is expressed in newtons 
For a trial design, set 


P EU “ P 4 - 12 N 


then 


R = 3.39 cm 

t . = 0.425 mm 

During stowage, the spiral material must elastically deform 
with a nominal strain of 


c 



(12) 


The allowable nominal packaging strain is a cr-^licated function of 
materials and detailed fabrication methods. We assume for prelimi- 
nary purposes an allowable value of 0.005 for graphite/polyimide. 
The trial design above gives 


t l _ 0. 425 

2R “ 2 x 33.9 


0.00627 


which is slightly too large. 

From Eqs. (9) and (10) we get 
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- = 4.14x10 (P 4 ) (P £U ) 


-x/6 


(13) 


For t /2R = 0.005, and P = 12 N. This gives 

V v 

P EU = 46 ’ 7 N 

R = 5.33 cm 

= 0.53 mm 

The strength test described earlier demonstrated that sufficient 
support is given to the longerons by the spirals if their cross 
section is half-round with a diameter of 2t^. For this case 


r 


A 


n 

2 


and the basic mass of an 8-meter-long batten is 87.6 grams. To this 
must be added masses to account for the finger hinges and end 
fittings . 
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Note: Each half batten molded on 7.11-cm dia 

mandrel from graphite/poly imide roving. 

Figure 2. Sketch of batten with dimensions showing a 
segment of lattice material (flattened) and 
a cross section of assembled batten. 
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Figure 4. Structural test model. 
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INTRODUCTION 


The Heliogyro concept was created in the mid-1960's by Astro 
Research Corporation (ARC) and MacNeal-Schwendler Corporation (MSC) 
under support of NASA Headquarters. References 1 and 2 summarize 
most of the results of that work. 

ARC has a subcontract with MSC to support JPL in applying the 
Heliogyro concept to the Halley Comet Solar Sailing Mission. 

William Ruff is the Technical Project Manager at JPL. The work at 
MSC is being led by Richard H. MacNeal, Program Manager. At ARC 
the work is being led by Karl Knapp, Program Manager. John M. 
Hedgepeth is ARC'S Senior Scientist. 

There are 12 blades in the Heliogyro design, and each blade is 
envisioned to be 8 meters in width and 7500 meters in length. The 
blades are expected to be composed primarily of a thin membrane 
constructed of material such as Kapton film with an aluminum reflec- 
tive coating on one side and an infrared emissive coating on the 
other. An overall view of the baseline Heliogyro is shown in 
Figure 1. 

This report presents the results of efforts made to eliminate 
the built-in wrinkles and pretensioned panels in the baseline design. 
The design was based on the data obtained from scale modeling per- 
formed at ARC and incorporates the results of previous analyses 
concerning wrinkles and membrane tensioning (References 3 and 4) . 


1 










c. 


ABSTRACT 


In the baseline Heliogyro design, wrinkles and pretensioned 
film panels control the dimensional changes in the film which ac- 
company thermal variations. However, wrinkles in the film affect 
both vehicle performance and dynamic stability. Thus, another 
method of tensioning and supporting the film panels was sought to 
allow for the expansion/contraction of the film material without 
the presence of wrinkles. 

A program of model making was initiated, coupled with a review 
of the previous analysis concerning wrinkles and sail panels. The 
experience gained in this program suggested a new panel-support 
design (see Figure 2). This design consists of biaxially-tens ioned 
film panels with parabolic edge members supplying the distributed 
constant edge loadings. 

The support system for each panel extends underneath the 
adjacent inboard and outboard panels (and vice versa) producing 
a high sail-to-vehicle area ratio. This interlaced design is 
biased towards the center of revolution to compensate for the 
radially varying centrifugal loading of each panel. With suitable 
variations of panel lengths, edge member tensions, and pivot point 
locations, this design will be applicable to ail radial stations 
of the 7500-meter Heliogyro blade assembly. 
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DISCUSSION 


The requirement for no wrinkles in the Heliogyro film panels 
has several far-reaching effects on the design of the attachment 
method for the film. This can be understood best by considering 
the reason for the presence of wrinkles in the baseline design. 

The panels had to be assembled with built-in sag and pretensioned 
chordwise strain to accomodate displacements resulting from thermal 
expansion and contraction of the film during a mission. The effects 
of not compensating for this strain are explained in Reference 4 
and consist of the unloading .. f the edge members during film con- 
traction and the production of chordwise wrinkling during expansion. 

Although the baseline design can accommodate the thermal strains 
of the film panels without significantly changing the stress dis- 
tribution in the blade assembly, there is an impairment in perform- 
ance of the Heliogyro from 5 to 10 percent due to the presence of 
the compensating slack in the film panels. The clack and accom- 
panying wrinkles also have undesirable dynamic characteristics 
because they may introduce periodic loading patterns near the 
natural pitching frequency of a blade assembly. The slack in the 
baseline design was necessary due to the difference in the expan- 
sion coefficients of the film membrane and the graphite edge mem- 
bers. For the projected thermal extremes, the graphite polyimide 
may be considered invariant, while the Kapton will change by ±0.75 
percent in the chordwise direction and ±0.69 percent in the spanwise 
direction. Thus, unless relative motion is allowed between the 
graphite polyimide edge members and the film, wrinkles are inevit- 
able. 

Another separate, though related, effect investigated in the 
modeling study was the problem of edge curl. If the material forms 
a closed cylinder at the free edge of a tensioned membrane, actual 
disintegration of the material may occur in the "solar oven" created 
by the curl. Incoming radiation is trapped within the dead space 
enclosed by the curl so that the temperature of the curl exceeds 
the melting point of the Kapton. If the edge of the film is rein- 
forced with additional material, the curl can be controlled as in 
the designs which incorporate tension members bonded directly to 
the Kapton. However, bonding a dissimilar material to the Kapton 
ciaates problems of differential expansion. Therefore, a design 
was sought which supports the panel and eliminates all conditions 


where curl and wrinkles might occur. If the material is cut away 
at a free edge in a parabolic scallop between supports, previous 
analysis has shown there will be no curl (Reference 4) . 

Thus, the requirements for the new design were an absence of 
wrinkles under all thermal conditions while maintaining acceptable 
stress levels in the film panels. Previous analysis was available 
describing wrinkles in tensioned membranes, but its application to 
Heliogyro panels had to be experimentally demonstrated. This is 
due to factors, such as a high aspect ratio, a varying body loading, 
and extreme temperature variations, whi~h make a straightforward 
analysis of the wrinkle properties of the Kapton difficult. The 
starting point of the modeling program verified the worst case 
indicated by the previous analysis (Reference 3) . This case showed 
that, whenever one of the principal mei/brane stresses vanished, 
wrinkling occurred parallel to the line of principal stress (see 
Appendix A, Figure A-l (a) ) . If the chordwise and spanwise stress 
are nonzero and positive, wrinkling should not occur. Practically, 
however, we find that if the ratio of one principal stress to 
another is large, say 50:1, wrinkling is possible. 

Using the values N x for chordwise stress and Ny for spanwise 
stress, one of the objectives of the test and modeling program was 
to investigate the value of N x /N y at the threshold of wrinkling. 

This number would be helpful in the Heliogyro design because it 
would help minimize the stress level in each panel which would, 
in turn, lengthen film life and reduce long-term creep. 
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MODELING 


A detailed account of the modeling process is found in Appen- 
dices A and B. A system of attaching the edge members to the film 
material was developed which allowed relative motion. Variations 
of woven, preformed, adjustable, and catenary-type edge members 
were made, resulting in an applicable system designed for solar 
sail requirements. (Typical models are shown in Appendix A, Figures 
A-l through A-3.) 

These preliminary models were made of 0.0001-inch mylar with 
an aluminum coating similar to the actual sail material. As the 
models became more representative of a feasible support system, 
a larger 1- by 2-meter panel was constructed out of 0.00025-inch 
aluminized mylar. This model had parabolic top and side edge mem- 
bers attached by flexible mylar loops bonded to the panel through 
which the polyimide edge members were threaded. 

The lower edge of the model was loaded by a Whiffletree 
arrangement which distributed a concentrated central load to the 
entire lower edge of the panel. By individually changing the loads 
on the edge catenary members and the bottom edge load, different 
values of N x , Ny, and N x /Ny were obtained. Table B-l in Appendix 
B lists the loadings which were studied. Because of the textured 
nature of the 0.00025-inch mylar and some creases which were inad- 
vertently created in the panel during assembly, all the photographs 
of the 1- by 2-meter panel show a partially wrinkled surface. 

Consequently, these photographs (shown in Appendices A and B) 
do not adequately reflect the degree of flatness which some loading 
points represent. A system of photographing the panel was developed 
which used backdrops, flash lighting, and a reflected grid back- 
ground (see Appendix B) . Even with this special effort to high- 
light the wrinkle patterns developed (or the lack of them), the 
film properties of this specific model always adversely affected 
the clarity of the photographs. Thus, in Figure B-12, Appendix B, 
there is no overall wrinkle pattern, which indicates a flat panel. 
However, the four distinct patterns of wrinkles shown in the photo- 
graph appeared throughout the testing and were inherent with the 
material used. Nevertheless, as much use was made of the 1- by 2- 
meter model in the time permitted. 

One particularly interesting pair of photographs, having the 
same edge member tensions, are shown in Appendix B, Figures B-l 
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and B-2 . The difference in the wrinkle patterns is due to the 
adjustment of the edge catenary, resulting in a slightly greater 
horizontal, or chordwise, stress (N x ) applied to the model. The 
approximate values of N x may be computed by the expression 

N = T/R 
x 

where T is the tension of the edge member, and R is the approximate 
curvature over the arc length. For this case, 

R = 10 meters 

T = 2.56 N 

N = 0.256 N A n 

x 

The effect of a 1-cm increase in the arc depression at this load 
could only account for, at most, a 20-percent increase in N x . Since 
it did not appear that this alone could be responsible for the 
presence of so many wrinkles or such a well-developed wrinkle pat- 
tern, another factor was assumed to be operative. Further experi- 
mentation and repetition of this test point was performed. Although 
the edge member of the catenary is flexible and conforms to the out- 
line of the catenary in the sail, after it emerges from the sail it 
must follow the tangential extension of the curve defined in the 
panel cutout. If it does not follow this extension, a constant dis- 
tributed load transmitted from the edge member into the panel will 
not be realized. Thus, once the C'.rve is cut into the panel and 
the edge member is properly aligned with the tangents at the end 
points, the only adjustment possible without changing the load 
distribution is the tension in the edge members. 

Another factor which affects this situation is the elastic 
properties of the panel itself. That is, as Ny increases, the 
length of the panel changes, much as the full-size panel will 
change due to thermal variations. For edge members of this design 
to work under different loads and dimensional variations, some 
allowance must be made for continuous alignment of the tension 
members to the panel curve. If the length o* a free edge member is 
short in comparison to the length threaded through the panel, a 
mechanical means of moving the edge member, relative to the panel, 
would have to be devised. However, if the point of support for the 
edge member is distant, relative to the panel, small variations in 
the location of the end of the panel will not affect the line of 
action of the edge member due to the small angles involved. 
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These factors explain the reasoning behind the long stringers and 
edge members discussed in the suggested alternative design covered 
in the next section. 

In regard to the determination of the threshold rat’o of N x /Ny 
for wrinkling to occur, this was found to be too complex to be ex- 
haustively examined in a test program of this scope. In general, 
if this ratio was greater than 0.1, the wrinkles present could be 
eliminated by adjustment of the edge members. For ratios less 
than 0.1, the presence of wrinkles was not only a function of N x /Ny, 
as the absolute load level of either load also seemed to be a fac- 
tor. As the loads became higher, the method of adjusting the 
position of the edge members and the precision of the layout of the 
edge members became more critical. 

In addition, the low aspect ratio of this model did not allow 
the end effects to distribute themselves over the width of the 
panel. This could have been partially alleviated by increasing 
the number of attachment points of the Whiffletree, but time did 
not permit this. Thus, in Figure B-15 (Appendix B) , although 
N x /Ny = 0.02, it is still difficult to tell if the wrinkle pattern 
developed is due to the overall loading, or just the higher stress 
concentrations in the regions of the lower attachment points. 
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SUGGESTED DESIGN ALTERNATIVE 


The results of the modeling indicated that, for a film support 
system to prevent wrinkles and still maintain a low level of stress, 
the positions and alignment of the supporting members had to be 
accurately placed. The alternate design presented by Astro Research 
satisfies the geometric and kinematic requirements for the support 
of the edge members (see Figure 2 and the enclosed sketch, SK1868) . 

The biaxially tensioned panel design, as seen in Figure 2, 
tensions the film such that no wrinkles will occur under any thermal 
variation encountered by the Heliogyro. Each panel is supported 
at each edge by parabolic tension members designed to impart the 
correct distributed load for that region of the panel. The use of 
interlocking parabolic end members allows for a minimum expansion 
gap between sail panels, as the ends of both panels are parallel 
and rectangular. 

The elements of the single-panel support system, shown in 
Figure 2, are: 

Item Desc r iption 

A* Collapsible lattice battens 

B* Outer edge members: tension varies from minimum 

at the tip to 650 N at the root 

C Inner panel edge members: tensioned at 42 N, 

1-mm in diameter, and constructed of polyimide 
graphite 

D* 2-pm film material 

E Polyimide graphite flat ribbon: connects film to 

parabolic tensioned members 

F Inboard panel end parabola: supports centrifugal 

load of panel 

G Outboard panel end parabola: maintains minimum 

tension along outboard end of each panel 

H Location of negator retractor for "G" : maintains 

constant tension for 0.75-m extension (see SK1868, 
Detail A) 

I Panel edge reinforcement 

♦Existing member in baseline design. 
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This arrangement of supporting members produces a minimum edge 
tension of 0.0125 N/m at each edge of the panel. Because of the 
slight chordwise contraction of the outer edge of the outboard end 
parabola, the stress in the center portion of the outer end of each 
panel will be reduced by approximately 50 percent. However, the 
stress ratio of N x /Ny will not be radically affected, as the span- 
wise stress, Ny, is also at a minimum here. 

Similarly, at the inboard edge of each panel, the spanwise 
tension is the greatest, varying from panel to panel depending on 
the radial station of the panel. In no case is it less than 
0.0125 N/m, and it is a maximum of 0.7 N/m at the tip. Since 
0.7 N/m t 0.0125 N/m is approximately 56, a shorter panel length is 
suggested for the more distant outboard panels. 

The expansion and contraction of each panel during temperature 
changes is taken up in large part by the negator spring which regu- 
lates the tension in the outboard parabola. The inboard end of 
each panel remains stationary relative to the edge members during 
expansion and contraction of the panel (except for chordwise motion) . 
As the panel expands, the looped edge members allow it to slide, 
relative to the long, edge catenaries. The chordwise changes in 
the panel are accomodated by the large-radius, small-angle devia- 
tion of the edge catenaries, which are tensioned by a rocker arm 
scissors mechanism connected to the main edge members (see SK1868, 
Detail A) . The displacements and loads for all temperature extremes 
have been calculated and appear in Appendix C. 
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SUMMARY AND CONCLUSIONS 


Because of the difficulty of treating the problem of wrinkles 
from a strictly analytic viewpoint, a series of scale-model film 
panels were r»dde. Various configurations and methods of attach- 
ments were attempted, and ultimately a larger 1- by 2-meter panel 
was constructed. The relationship between stress levels, precision 
of application, and precision of loading was more apparent in the 
larger model. Once a satisfactory method of supporting the sail 
panel had been obtained, the problem of thermally-induced dimen- 
sional variations was examined and a promising solution has been 
outlined. 
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Figure B-6. 

Loading #6, N x /N y = 0.047 
(some wrinkles appear at 
lower corners) 


Figure B-5. 

Loading #5, N x /N y =0.10 
(stringers placed in 
front of panel) 
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Figure B-7. 

Loading #7, N x /N y = 0.040 
(more wrinkles appear 
at lower edge) 


Figure B-8. 

Loading #8, N x /N y =0.17 
(first attempt 
at reflective grid) 
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Figure B-ll. 

Loading #31, N x /Ny = 0.04 
(same as Load #7 with 
grid added) 
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1. INTRODUCTION 


A significant weight penalty (of the order of 300 kg) is paid against 
the micrometeoroid hazard in th . baseline design. Estimates of the amount 
of weight required vary by a factor of four or more, because of uncertainties 
regarding the damage caused by meteoroids. There are uncertainties rega:ding 
the physical distribution (density, velocity and flux) of meteoroids, but 
even larger uncertainties exist regarding the damage caused to an edge 
tendon by a meteoroid with given size, density and velocity. It is proposed 
that, in order to reduce the range of uncertainty, a research effort be 
undertaken which concentrates on the damage to one-dimensional structures 
similar to those proposed for the edge tendons of the Heliogyro, using 
current estimates of the physical distribution of meteoroids. 

The study should employ both analytical and experimental methods. 
Experiments are necessary because the available experimental data was 
obtained at unrealistically low velocities, and because no data at all is 
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available for the critical case of low angle (grazing) impact, which is the 
major source of uncertainty in the meteoroid hazard. 

The results of the study will be useful in the future for the design of 
any large space structure. The most efficient way to carry a small load for 
a long distance is to use a framework of extremely slender tension and 
compression members, rather than panels. In the case of large structures 
for space, the member sizes tend to be in a range where meteoroid damage is 
an important consideration. Examples: Heliogyro edge tendons and battens, 

square-sail ties. Astro masts. 

2. REVIEW OF CURRENT STATUS 


Although a literature search has not been made, it is the opinion of 
qualified personnel at JPL that no work on micrometeoroid damage to fila- 
mentary structures has been published between Ref. 1 (1967) and the present 
time. 

Reference 1 defines a procedure for calculating the probability of 
failure of a filamentary structure due to micrometeoroids. The procedure is, 
in part, based on the experimental data that was available at the time and, 
in part, on assumptions regarding physical behaviour which are unsupported 
by experimental data. The critical assumptions are in regard to the extra- 
polation of hole size for normal particle incidence from experimental 
velocities (2-8 km/sec) to meteoroid velocities (-30 km/ sec) , and in regard 
to the damage caused by grazing incidence. 


The procedures of Ref. 1 were used at MSC and at Astro in January 1977 
to estimate meteoroid damage to the edge tendons of the Heliogyro. The 
estimates indicated that the probability of failure was less than l/10th 
percent for a trifilar edge tendon with 265 kg total weight. 
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Later, Bamford at JPL discovered numerical errors in the procedure of 
Ref. 1 for the extrapolation of hole size for normal incidence, and sub- 
stituted an extrapolation proposed in Ref. 2, which gives much larger holes. 

He retained the assumptions of Ref. 1 regarding grazing incidence and used 
a current estimate of meteoroid flux distribution which is more severe 
than that in Ref. 1. With these modifications to the procedure, he found 
that the trifilar design has nearly a 100% probability of failure. He then 
proceeded to develop an edge tendon design with five .001 in. tapes that 
weighed about twice as much as the trifilar design, and which he found 
(Ref. 3) to be barely acceptable (4% probability of failure). 

Recently, MacNeal (Ref. 4) has proposed another method for extrapolating 
the hole size for normal incidence which gives much smaller holes than the 
method used by Bamford, but which also fits the available experimental data. 

In any case, it may be concluded that we really don't know how to extra- 
polate .he hold size for normal incidence-, but that it makes a big difference. 

The situation with respect to grazing incidence is even worse. The 
basic assumption of Ref. 1 was simply to take the hole size for normal inci- 
dence and divide it by the sine of the incidence angle, a procedure which, 
it must be said, is quite reasonable for particles that are very large 
compared to the thickness of the sheet. Calculation showed, however, that 
a mathematical singularity (i.e., a 100% probability of failure) exists for 
an edge-on hit. The singularity was removed by the semi-rational assumption 
that a particle striking the tape at a small incidence angle will break the 
tape only if it can also break a solid round wire with the same cross- 
sectional area. Even so, calculations based on the theory show that most 
of the damage is done by very small particles at grazing incidence angles. 
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Justification for the grazing incidence theory of Ref. 1 is slender 
indeed. Experimental and analytical evidence presented in Ref. 2, show, 
fairly conclusively, that meteoroids traveling at 30 km/sec will be vaporized 
when they strike a target, even a very thin one. If this is so, then it 
couid be reasoned that the back end of a high-velocity particle striking a 
target at low incidence will be vaporized before it reaches the target and 
that the resulting damage to the tape will be surface scorching which 
spreads out from the initial point of contact, rather than an elongated 
hole. In this case, failure of the tape could be computed from the im- 
pulsive load caused by stopping the normal momentum component of the particle. 
Such calculations have not been made. 

Again, we don't know much about grazing impact but it makes a big 
difference. For example, consider the choice of flat tapes vs. curved tapes. 
If meteoroids simply continue in a straight line and make holes just big 
enough to pass through, then curved tapes are better than flat tapes 
because the problem of extremely elongated holes at low incidence is 
avoided. On the other hand, if the vaporization theory is correct, chen 
curved tapes might be worse than flat tapes because the gas released by 
the impact will cause higher pressures on the tape. In summary, the 
current status is that we lack a reliable procedure for calculating 
meteoroid damage. One is needed because, by current (unreliable) estimates, 
the weight required to counter the meteoroid threat is significant. 

3. PROPOSED RESEARCH PROGRAM 
3. 1 Experimental Program 

Facilities are available which can accelerate small particles to the 
meteoroid range (30 km/sec) . It is proposed that such facilities be used 
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with targets similar to those proposed for the Heliogyro edge tendons (i.e 


thin graphite-polyimide tapes). Many firings will be required to cover the 


ranges of the relevant parameters. The following parameter ranges are 


recommended 


1. Particle velocity (10 km/sec to maximum available velocity) 


Particle density (two or more values including one that is near 
0.5 gm/cm . This might be achieved with hollow spheres.) 


Target thickness (.001 inch). Target thickness can be varied 
instead of particle mass, if particle mass and velocity parameter 
ranges cannot be met. 


6. Target width (1 cm to 3 cm) 


7. Target length (10 cm or more) 


8. Target materials (graphite polyimide tape, aluminum tape) 


9. Tension in target material (zero, 1/3 ultimate stress, 2/3 
uicimate stress) 


10. Target chordwise curvature (included angle = 0 


Combinations of parameter values should be carefully selected to maximize 


the value of the knowledge gained for a given cost. This will require that 


some analytical work precede the experiments 


A Second set of experiments is recommended to validate proposed 


designs of tendon assemblies, by observing damage when a member is partially 


or completely cut 


It may also be necessary to study impulsive failures of individual tapes 
with a separate facility, where the impact is simulated by explosive charges 


or other means 


K, 
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3. 2 Theoretical Program 

A theoretical program is needed to provide empirical formulas that can 
extend test data to new cases and which can be used to design structures. 
Research in the following areas is suggested: 

a. Nature of the impact. Size and shape of the hole. Phase (solid, 
gaseous) and angular distribution of the debris. Surface damage 
when there is no hole. Research in this area involves physics, 
fluid dynamics and thermodynamics. 

b. Failure of tapes due to impulsive loads. Research in this area 
involves analysis based on structural dynamics principles. 

c. Design of redundant systems. How many tapes. How wide. How 
far between load transfer points. Research in this area involves 
probability theory and detailed stress analysis. 
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PROPOSED MODIFICATIONS OF NASTRAN TO IMPROVE HELIOGYRO ANALYSIS 


1. INTRODUCTION 


The NASTRAN computer program has been used in the Solar Sail design 
project to perform dynamic analysis of the Heliogyro. The dynamic charac- 
teristics of the Heliogyro are similar to those of a conventional heli- 
copter rotor, but the differences are sufficiently great that computer 
programs designed for conventional rotor analysis (such as the SADSAM 
program, Ref. 3) cannot be used without extensive modification. Although 
the NASTRAN computer program has no specific capability for rotor analysis, 
it is, in general, extremely versatile and easy to modify and has the 
required capacity. The decision to use NASTRAN for analysis of the Helio- 
gyro has been justified by the results that were achieved, but exper- 
ience in that effort also indicates the need to improve some aspects of 
NASTRAN to provide a more usable and less costly analysis procedure. 

This report outlines the tasks which are considered necessary to carry 
out the required analysis types (dynamic stability analysis and trans- 
ient response analysis) in a design environment. 


2. IMPROVEMENTS TO THE HESSENBERG METHOD OF COMPLEX EIGENVALUE EXTRAC- 


2.1 Background 


NASTRAN includes three methods of complex eigenvalue extraction, 
but the only efficient method for general use is a transformation method 
called by the keyword "HESS." In this method, the second order equations 




of dynamic motion are converted to first order form and are then trans- 
formed by Hessenberg reduction to upper Hessenberg form. The QR-method 
is used to compute the eigenvalues, and inverse iteration is used to 
compute the eigenvectors. 

The present version suffers from the following limitations: 

1. All coefficients must remain in memory. In prior Heliogyro 
work, this led to a limitation of about 40 dynamic degrees of 
freedom. 

2. The logic for limiting the number of eigenvectors to be com- 
puted is deficient, requiring that all be computed if any are 
needed. 

3. The solution is sensitive to matrix scaling. The Heliogyro 
has low natural frequencies (.005 cycles/sec). It was found 
necessary to scale the equations of motion so that the lowest 
natural frequency was near 1 cycle/unit time in order to obtain 
reliable solutions. 

4. The mass matrix is inverted, requiring it to be nonsingular. 

The solution becomes less reliable if it is nearly singular. 
While this limitation can always be avoided, the resulting 
modeling techniques become tedious and nnecessarily compli- 
cated. 

2 . 2 Technical Approach 

The size limitation will be expanded by rewriting the algorithm to 
release unneeded space in memory. The user will have the ability to 
define the region in eigenvalue space where eigenvectors are to be deter- 
mined. Matrix scaling will be done inside the solution process, making 
it transparent to the user. Two options for allowing singular mass 
matrices will be investigated. 
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2.2.1 Relaxing Size Limitations 


The upper Hessenberg matrix has nonzero terms on and above the 
diagonal, and on one subdiagonal. It is proposed that this matrix be 
transposed, and given another Hessenberg reduction, so that it becomes 
both upper and lower Hessenberg in form, i.e., tridiagonal. The eigen- 
values and eigenvectors of tridiagonal matrices with several thousand 
degrees of freedom can be solved with all coefficients in memory. Since 
this is far beyond the requirements' of the Heliogyro analysis, the only 
operation that requires further consideration is the Hessenberg reduction 
itself. Wilkinson (Ref. 1) has stated that the second Hessenberg reduc- 
tion has some theoretical shortcomings with respect to numerical stab- 
ility, but that in practice these become evident on only a small range of 
pathological problems. 

The basic Hessenberg reduction is defined by the equations 

[A] [N] = [N] [H] (1) 

(See Ref. 1, p. 355-412.) [A] is the matrix whose eigensolution is to 

be found, [N] is a unit lower triangular matrix, and [H] is the upper 
Hessenberg form of [A] . 

Reference 1 describes a technique for performing the reduction in 

memory by overwriting terms of [N] and [H] on [A] as the terms are com- 

2 

puted. This results in a storage requirement of 2(N A + N^) words for 
single precision, complex matrices where is the number of columns in 
[A]. Since the present method requires 6N“ ♦ 8N A words, considerable 
improvement can be made. This would raise the present limit of about 40 
degrees of freedom to about 70 degrees of freedom for computers similar 
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in capacity to the Univac 1108 used in previous Heliogyro work. Seventy 
degrees of freedom (modal coordinates) is regarded as adequate for Helio- 
gyro work. Higher capacity is available on larger computers. 

The present QR iteration for eigenvalues requires storage of [Q] , 

a full unitary matrix, [R] , an upper triangular matrix, and [A], which 

2 

is almost triangular in form. This results in 4(N^ ♦ N^) words for single 

precision, complex matrices, which would become the new limiting factor. 

However, as mentioned above, another stage of Hessenberg reduction will 

T 

be performed on [A] , so that both it and the [Q] and [R] matrices will 

retain their tridiagonal form throughout the iteration. Thus, the 

2 

required storage will be that for the Hessenberg reduction, 2(N^ + N^) . 

At some future date, consideration will be given to the addition of 
spill logic to HESS. This will remove the limitation on problem size 
imposed by available core memory, but the practical limit imposed by 
cost considerations is at present not much larger than seventy degrees 
of freedom. 

2.2.2 Restricted Eigenvector Calculation 

The user will control the number of eigenvectors calculated and 
output by the method described in remark 8 on the modified EIGC Bulk 
Data card shown in Figure 1. 

2.2.3 Scaling of Coefficient Matrix 

The [A] matrix will be equilibrated to avoid scaling problems, 
using the method described in Ref. 1, pp. 356-357. The basic equation is 





w a2 

a b2 

l 'o2 


o.3 

5.6 

5.6 

1.5 



Contents 

Set Identification number (unique. Integer > 0). 

Method of complex eigenvalue extraction, one of the BCD values "INV," “OET," or 
•HESS.* 

INY - Inverse power method 
DET - Determinant method* 

HESS - Upper Hessenberg method 

Method for normalizing eigenvectors, one of the BCD values "MAX" or "POINT." 

MAX - Normalize to a unit value for the real part and a zero value for the 
Imaginary part the component having the largest magnitude. 

P0INT - Normalize to a unit value for the real part and a zero value for the 
Imaginary part the exponent defined in fields 5 and 6 - defaults to 
"MAX" If the magnitude of the defined component Is zero. PCINT Is not 
available for modal formulations. 

(Continued) 

Figure 1. Modified EIGC Bulk Data Card. 
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NASTRAN DATA DECK 


EIGC (Cont.) 


G . 

C 

E 

( °«J‘ V 
l J 


Grid or scalar poi.it identification number (required If and only if N0RM»P0 1 NT ) 
Integer >0). 

Component number (required if and only if N0RM-P0INT and G is a geometric grid point) 
(0 < Integer < 6). 

Convergence criterion (optional) Real > 0.0) 

Two complex points defining a line In the complex plane (Real) 

Width of region In complex plane (Real > 0.0) 

Estimated number of roots in each region (Integer > 0) 

Desired number of roots in each region (Default is 3N g ^ ) (Integer > 0) 


Remarks : 1. The preferred method is HESS, provided that sufficient main storage is available. 

Insufficient storage for HESS will cause the program to switch to INV. 

2. Each continuation card defines a rectangular search region for METH0D*INV or DET. 
Any number of regions may be used and they may overlap. Roots in overlapping 
regions will not be extracted more than once. 

3. Complex eigenvalue extraction data sets must be selected in the Case Control Deck 
(CMETH0D»SID) to be used by NASTRAN. 


4. The units of a, u are -aUio.'.s per unit time. 

5. At least one continuation card is required. 

6. For the determinant method with no damping matrix, complex conjugates of the roots 
found are not printed. 

7. See Section 10.4.4.5 of the Theoretical Manual for a discussion of convergence 
criteria. 


8. For the Upper Hessenberg Method, N „ controls the number of vectors 

computed. All vectors for roots in the area defined by the (a, w) 
pairs will be completed and output, if these values are input, 
regardless of the value of N^. 


9. The required working storage for the Upper Hessenberg Method Is given by 

U * 6N 2 ♦ 8N 

where N is the order of the stiffness matrix if thrre is no damping, and twice the 
order of the stiffness matrix if damping is present. 

10. If Method * MESS- the mass matrix must be nonsingular. 


Figure 1. Modified EIGC Bulk Data Card (Cont.) 
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W„ew- f D J t A] o i df D g* 1 C2) 

where f D J is a diagonal matrix. 

2.2.4 Allowance for Singular Mass Matrices 

Two techniques will be investigated. The easiest to implement in- 
volves shifting the eigenvalues of the problem. Instead of inverting 
the mass matrix [M] , the shifted matrix [p^M ♦ p s B + K] is inverted. 

This allows singularities in [M] tobe repressed by terms in the other 
matrices. This technique has proven effective in the modified Givens 
method in MSC/NASTRAN. 

An elegant but more difficult method identifies singular degrees of 
freedom in [M] , and reduces the size of the problem in a manner analogous 
to static condensation, (see Ref. 2). Although this method may cost 
more to implement, the reduction of problem size operates favorably on 
the cubic cost curve typical of transformation methods, by suppressing 
the uninteresting infinite frequency modes of the system. 

3. IMPROVEMENTS TO THE TRANSIENT RESPONSE CALCULATION MODULES 
3 . 1 Background 

The transient response modules are the largest consumers of computa- 
tion time for HeliogyTO analysis, and require the following new capabi- 
lities in addition to general efficiency enhancements: 

1. Improved spill logic: the ability to solve problems too large 

to fit into memory by placing some portions of the calculations 
in secondary storage, if necessary. 
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2. Rotating to nonrotating coordinate transformations ("Resolver"): 
The method of Lagrange multipliers has proven effective in 
prototype work at MSC and will be added to the TRD module. 

3. Nonlinear, timevarying coefficients: certain structural and 

photodynamic terms must be updated at every time step. The 
only practical method to accomplish this is by placing the 
equations that generate these terms inside the transient 
response calculation. 

4. Data Reduction: the Floquet method of system identification 

will be implemented to recover the frequency and damping of 
system modes from the time histories of transient response. 

Numerical integration procedures must be carefully designed to 
treat problems with few dynamic degrees of freedom inexpensively while, 
at the same time, imposing no absolute limit on the solution of large 
problems . 

The computation of transient response may be divided into two 
phases: 

a. Preparation phase: assembly and reduction of the equation of 

motion; 

b. Numerical integration phase. 


The key to a low-cost solution is to do as much work as possible in 
the Preparation Phase in order to limit the amount of data accessed at 
each integration step. Tne reason is that, if the amount of data cannot 
fit within the available high speed memory, the cost of accessing it 
will be the dominant element in the total cost. Thus, the problem size 
for low-cost solutions has a practical limit which depends on the size 
of the available high speed memory. Careful planning of the data pro- 
cessing procedures is required to ensure that this critical size will be 
as large as possible. 


- 8 - 


The NASTRAN procedures for numerical integration need to be re- 
viewed and upgraded. NASTRAN tends to emphasize the solution of pro- 
blems with very many degrees of freedom rather than low cost solutions 
of moderate uize problems. Our SADSAM program has the opposite em- 
phasis. It produces very low cost rotor analysis (Ref. 3) but problem 
size is limited by fixed dimension statements. What is really wanted is 
both capabilities, however incompatible they may appear to be. 

With regard to integration algorithms, NASTRAN currently provides a 
choice of two -- a version of the Newmark Beta method for general use 
and a semi-analytic recursion method for uncoupled modal equations. We 
use the same version of the Newmark Beta method for rotor analysis in 
SADSAM with considerable success. 

With regard to the treatment of nonlinear effects, it may be as- 
sumed that each nonlinear force or moment is a function of the motions 
at particular points. The matrix which relates these motions to the 
reduced dynamic freedoms will be core held, in packed form, during 
numerical integration. The user may reduce its size either by reducing 
the number of reduced dynamic freedoms, or by limiting the number of 
points to which nonlinear forces and moments are applied. The F0RTRAN 
code which computes nonlinear forces and moments, and the associated 
tabulated empirical data, will also be core held during numerical inte- 
gration; so will the integration matrices which are derived from the 
linear mass, damping and stiffness coefficients for the reduced dynamic 
freedoms. If all of the data cannot fit in the high speed memory, some 
of it will be kept in disk storage and transferred to the CPU at each 
time step. The solution vector, plus other data needed for post-processing, 
will be transferred to disk storage at each time step. 
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3.2 Technical Approach 


3.2.1 Design Goals for the Transient Analysis Modules 

A "standard"-size Heliogyro model is defined here to consist of one 
blade with a fully or partially fixed hub. The transient response 
module will be modified to solve this model without spill. A "large"- 
size Heliogvro model is defined as one with six blades, payload, and 
control system modeling. The module will be designed to solve this 
problem routinely, using spill logic if necessary. Upper limits on 
problem size are tabulated below. 


Symbol 

N 


Meaning 

noA. 

Standard 

Large 

Grid points x 6 d.o.f. each 

300 

2000 

Dynamic variables (not eliminated 
by constraints or static conden- 
sation) 

150 

1000 

Forcing points 

75 

500 

Modal variables 

15 

80 

Time steps 

1000 

1000 


3.2.2 Expansion of Present Modules 
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Time-dependent loads are presently generated and reduced to modal 
generalized forces {P^} in the TRLG modi le. This module is adequate for 
HeliogyTO anal) sis. So is the numerical integration performed in the 
TRD1 module. 


A modified version of the TRD1 module will be provided for Helio 
gyro analysis. The basic integration algorithm will be 

[DHu.^I = {P 4 ^,} ♦ {N ;i1 } ♦ [C] {u. _ , > ♦ [E] {u. } ♦ (N . , } (3) 





The only term not used in the present algorithm is {N^ + ^}. [C], [0] and 

[E] are constant matrices for a given time step site. The {u^} vector 
represents the modal variables {u^} at the time step i. {rr + ^} is 
a nonlinear term (i.e. f displacement or velocity dependent) which is 
available now in NASTRAN but has -estricted utility. 

The new term {N^ + ^} represents nonlinear, time-dependent forces 
applied directly to grid points. Two types of nonlinear forces that are 
useful for Heliogvro analysis are described in Refs. 4 and 5. The non- 
linear forces are generated by the following equations, performed at 
every time step: 


(u 


d, i+1 ' 


lp n4 } = fft 
‘ d,i*l J v i* 


{N i+1 > = [4> 



(4) 

'“d.i 1 - {u d,i.l» 

(S) 

] T {p ni } 
dh J 1 d,ifl ; 

(6) 


The [0^] matrix is a compressed matrix of eigenvectors whose coeff- 
icients for u, variables without nonlinear forces are set to zero. It 
a 

is stored in packed form in single precision. The {u^} variables are 
physical motions at the nonlinear stations, represents the 

resulting nonlinear physical forces, and {N^^} is the vector of gen- 
eralized modal nonlinear forces. 


Memory Requirements 



Equation Number 

Terms 

Sizes 

(3) 

rci. [D] , [E] 

3N? 

n 

(3) 

{u i } ’ { Vi } - {F W* {N i } * { *i } 

5N, 

n 
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Equation Number 


Terms 


Sizes 


(4) 

(5) 

( 6 ) 


{u d,i*l K [ *dh ] 


{P ni } 

1 d,i+l' 


{ *i } 


N, 



N 


N 


P 

h 


where is the dimension of {u^} and is the number of forcing points, 

2 ( *1 

Total storage requirement: 3N^ ♦ 6N^ ♦ N ♦ (N^ x N ) . 

CPC IBM/Univac 

Single precision words for standard problem 1,200 1,602 

Single precision words for large problem 60,000 80,000 

The high-speed memory requirement for large problems is within the 
range of possibility on computers available today. For example, on the 
Univac computer used on past Heliogyro analysis, there are 42,000 words 
of core used for code, and 85,000 words of open core available for data 
storage with the TRD1 module. After the data center which operates the 
computer completes a scheduled hardware upgrade, there will be 127,000 
words of open core available for data storage. Spill logic will be 
provided for the module, but it appears that spill will not be required 
for Heliogyro analysis. 

3.3.3 Resolver Capability 

The equations of motion of the blades are most conveniently ex- 
pressed in rotating coordinates, while those of the control system must 
relate to a nonrotating reference frame. The time-varying transforma- 
tion will be modeled by the technique of Lagrange multipliers (see Ref. 

(*)Single precision on all computers. Other variables are double precis 
on IBM and Univac. 


At 

* 


i 



6, Section 4.6). All degrees of freedom in the nonrotating coordinate 
system will be sequenced to be the last n variables. The time-dependent 
terms will be added to 

3.3.4 Floquet System Identification 


The HeliogyTO modes have low natural damping, so that determination 
of system stability by inspection of transient analysis requires inspec- 
tion of hundreds of rotor revolutions before asymptotic behavior is 
achieved. Methods to determine closed-form frequency and stability 
parameters for the system based on the Floquet hypothesis have been 
described in the literature (see Ref. 7). The basic integration al- 
gorithm (Equations 3 through 6) is modified as follows: 


After the system has reached an interesting state, a time of 

initial obsex-vation t is established, based on user input 

(i.e., after y rotor revolutions, or after the variation in 

response between successive cycles is less than z percent) . 

Equation 4 is solved to find the state at time t . named 

o 


[D]{u i+2 > = {P} * [C]{u i+1 > «■ [E](u i ) 


( 7 ) 


[P] contains the time-dependent and nonlinear forces. 


A matrix of perturbed initial states 1 *- s determined. 


l . e. , 


t’hh'W 1 


U. _ u. - u. , 

i*2 i*2 i*2 


[IC] 


( 8 ) 


.An obvious choice for the perturbed initial condition matrix, [IC], 
is a diagonal matrix of small terms. Provisions will be made for 




of 
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user-selection of this initial condition and several others. The [C^] 
matrix is then integrated using Equation 7 above and the same startup 
techniques as the present algorithm. At the t n time step, Equation 7 
is of the form 

[D][$(t n> t o n = {P} ♦ [C][$(t nl .t 0 )] ♦ [E][$(t n2 ,t o )l (9) 

This equation is evaluated at every time step over one rotor cycle. 

Note that only three additional N. x matrices need be stored at one 
time. 

Reference 7 describes the method of computing the frequency and 
stability parameters of the system from the state transition matrix 
after one period. This matrix can be computed from 

[0(t o ♦ T,t 0 )] = [$(t Q * T,t o ) -1C] [IC]" 1 (10) 

where T is one rotor period. Small changes to utility routines will be 
needed to compute logarithms of complex diagonal matrices. The calcu- 
lations will be performed in a new module. 

The use of a reduced set of initial conditions for approximate 
answers is a subject of current research. For example, 11/2 initial 
vectors rich in the lowest modes would be expected to produce response 
dominated by the lowest modes. A least-squares fit can be used to find 
the reduced [d> (r q ♦ T,t Q )] matrix. As the cost of this operation 
is cubic with N^, cost and memory site reduction can be quite dramatic, 
and may even improve the accuracy of determination of the low-frequency 
modes. Provision will be made to use a variety of initial conditions 
tnd smoothing techniques. 


- 14 - 



t 




Since this is a research topic, an experimental version of NASTRAN 
will be used as a test bed for developing a production tool. Low cost 
but inconvenient input and output formats will be used until the prac- 
ticality of the method has been demonstrated. 

3.3.5 Summary of Enhancements for Module TRD1X 

New Input Data Blocks 

Eigenvector transformation 

UC0NTR0L General data block, input by user on DTI Bulk Data cards 
New Output Data Blocks 

[4>(T)] State transition matrix at end of one rotor cycle 
(optional) 

[<J(t)] State transition matrices at every time step over one 
rotor cycle, stored in appended form (optional) 

EXP0UT General output data block (optional). 

Method 


1. General input and output data blocks are used for data during 
module testing, and for unconventional features. Input for system 
identification research, such as specification of initial condi- 
tions, will be done on DTI, UC0NTR0L Bulk Data cards. Special 
debug output, or output requirements not anticipated, will pass out 
through the EXP0UT data block. These blocks allow module changes 
without requiring MPL updates. 

2. Inspect the new nonlinear data to determine unloaded rows 
[$dhl • Eliminate these rows and convert to single precision i r 
necessary. 
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3. Compute the new nonlinear forces and extend the integration 
algorithm to accommodate them. Upgrade the integration routine 
to store all needed data in core. 

4. Improve user control of output. Allow integration over a timf.- 
span without any output. 

5. Install code for resolver. User inputs will be on DTI Bulk 
Data cards. 

6. Install code for Floquet system identification. 

4. DMAP SEQUENCES FOR SOLUTION OF STRUCTURES IN ROTATING COORDINATE 

SYSTEMS 

4 . 1 Background 

Past Heliogyro analysis in NASTRAN was done with ad hoc DMAP alters 
to the NASTRAN rigid formats. This was effective during exploratory 
studies, when the significance of many second-order effects had to be 
determined. It resulted in twenty rigid format ALTER packages, many of 
which have only subtle differences. This exploratory wov : ; has shown 
which effects must be included in the analysis. This will allow con- 
solidating the technology into the general-purpose DMAP sequences listed 
below. 

A new capability has been developed in MSC/N'ASTRAN since the prior 
work was completed. It is a small strain, large deflection statics 
capability which accounts for geometric nonlinearity, including higher 
order terms not included in the present differential stiffness capa- 
bility. Use of this technology and the DMAP sequences listed below will 
decrease the labor and calendar time needed to assemble Heliogyro models, 
and will increase the portability of this technology to JPL. 


( 
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4.2 DMA P Seque nce to Compute Lin earized Structural Matrices from 

Geometric Nonlinear Analysis 

User Input - Initial geometry, rotation speed, sun load intensity, 
structural description, blade pitch angles. 

Output - Structural matrices for perturbations from the equilibr. 
position, bulk data cards or their equivalent defining the structure at 
its equilibrium position, the Coriolis force damping matrix, "tennis 
racket" effects (moments due to the angle between the rotor plane and 
the principal axes of inertia), and photodynamic influence coefficients. 

4 . 3 DMAP Sequence to Compute System Stability 

User Input - The matrices from the DMAP sequence described in 
Section 4.2, the hub constraints for hub-fixed modes, collective modes, 
and cyclic modes, plus models of damping devices and control systems. 

Output - Eigensolutions , stability margins, energy absorbed per 
mode by damping devices. Both modal and direct solution techniques will 
be available. The Lagrange multiplier technique will be used to imple- 
ment the complex constraint equations needed for cyclic modes. 

4 . 4 DMAP Sequence for Transient Analysis 

User Input - Linearized matrices from the DMAP sequence described 
in Section 4.2, control perturbations, nonlinear photodynamics , non- 
linear Coriolis effects and nonlinear solar illumination pressure forces. 

Output - Time histories of motion and internal forces, options for 
roots of perturbation solutions to the nonlinear transient solution 
using the Floquet system identification technique. 
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ESTIMATE OF THE REQUIRED EFFORT 


The following table presents an estimate of the manhours required 
to perform the tasks described above. It assumes that the work will be 
done by senior engineers and programmers who are experienced in NASTRAN 
development. Computer time is not estimated but experience with this 
type of development shows that computer cost is about one-half of the 
manhour cost. 


Task Manhours 

1. Rewrite code for the Hessenberg method of 550 

Complex Eigenvalue extraction 

2. Modify the transient response module 360 

3. Install the resolver capability 210 

4. Install Floquet System Identification 270 

5. Write DMAP sequtnces for Heiio^yro Analysis 410 

Total 1,800 


REFERENCES 


Wilkinson, J. H. , The Algebraic Eigenvalue Problem , Oxford Univ. 
P-ess, 1965. 

Dell, A.M. , Weil, R.L., Thompson, G.L., "Roots of Matrix Pencils: 

The Generalized Eigenvalue Problem," Communications of the ACM, 

Vol . 14, No. 2, pp. 114-117, Feb. 1971. 

SADSAM User's Manual , MacNeal-Schwendler Corp. Report MSR-10, 

Dec. 1970. 

Bellinger, E.D., ,r Modification of NASTRAN for the Nonlinear 
Coriolis Force Calculations of a Rotating Blade," MacNeal-Schwendler 
Corp. Memo EDB-2, Project EC-404, Heliogyro Preliminary Design. 

Bellinger, E.D., "Modification of NASTRAN for Nonlinear Solar 
Illumination Pressure Force Calculations on a Heliogyro Blade," 
MacNeal-Schwendler Corp. Memo EDB-3, Project EC-404, Heliogyro 
Preliminary Design. 

MacNeal, R.H., "The Dynamics of Rotating Elastic -Bodies," MacNeal- 
Schwendler Corp. Report MSR-36, August 1973. 

Gockel, M.A., "Practical Solution of Linear Equations with Periodic 
Coefficients," AHS Journal, Jan. 1972. 





— 



MS404-3 


CONCEPTUAL DESIGN OF A 
FLIGHT TEST MODEL OF THE HELIOGYRO 


Prepared by 


Richard H. MacNeal 


28 October 1977 


THE MACNEAL- SCHWENDLER CORPORATION 
7442 North Figueroa Street 
Los Angeles, California 90041 


-*■ 



» t 



TABLE OF CONTENTS 


Section 

1.0 

2.0 

3.0 

4.0 

5.0 

6.0 

Table 1 

Table 2 

Table 3 
Figure 1 
Figure 2 
Figure 3 
Figure 4 
Appendix 


INTRODUCTION 

DESIGN CRITERIA FOR A FLIGHT TEST PROGRAM . . 
DESIGN CONSTRAINTS AND SELECTION OF PARAMETERS 
DESCRIPTION OF THE 4.1 METER FLIGHT TEST MODEL 

PARAMETRIC VARIATION OF BLADE CHORD 

CONCLUDING REMARKS. . . .' 


Comparison of Parameters for the Full-Scale Halley Vehicle 
and the 4.1 Meter Flight Test Model 

Comparative Mass Distributions for the Full-Scale Halley 
Vehicle and the 4.1 Meter Flight Test Model 

Detailed Mass Distributions 

Four Meter Heliogyro, Stowed Configuration 

Four Meter Heliogyro Blade, Deployed Configuration. . . . 

Performance Comparison 

Sail Module Mass vs. Deployment Reel Length 

Design M ethod 




MS 404 -3 



1 

! 

i 

I 

i 

\ 

I 



Conceptual Design of a Flight Test Model of the Heliogyi • 

1.0 INTRODUCTION 

As part of the design review of the Heliogyro Solar Sail proposed for 
the Halley Rendezvous mission, a need was identified for a flight test 

program. This report is addressed to the conceptual design of a vehicle 

capable of providing the required information. The vehicle will also be 
able to perform useful planetary missions. 

2.0 DESIGN CRITERIA FOR A FLIGHT TEST PROGRAM 

Although the length of Heliogyro blades exceeds the size of any 
ground-based facility, the mechanical components involved in deployment 
are of comparatively modest dimension J8 m) , so that deployment tests 
can be performed in the laboratory. Special care in removing gravita- 
tional effects will not be required, because the centrifugal forces 
during deployment are of the order of 1/2 g. Thus, ground-based deploy- 
ment tests should be considered to be adequate and reliable. 

The only major area in which ground-based tests are not feasible is 
the area of dynamics and control, which involves complex interactions 
between photon pressure and centrifugal force. The forces due to photon 
pressure acting on a very thin film are of the order of .001 g so that a 
vacuum chamber would be required. The largest scale model that could fit 
in a vacuum chamber is of the order of 1/1000 scale, which is too small to 
provide useful information. Tims, only a flight test will be able to 
verify theoretical dynamic calculations. 
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In order to provide useful dynamic information, the flight test 
configuration and the environmental conditions should be reasonably 
similar to those of the full-scale vehicle. For this reason the follow- 
ing should be avoided: 

a. Low altitude flights where atmospheric drag is appreciable 
compared to photon pressure. 

b. A vehicle with fewer than four blades. (This is the smallest 
number that can provide dynamic similitude with the full-scale 
vehicle.) 

c. Very small vehicles, which will of necessity violate dynamic 
similarity, for example, by having a central hub which is 
much heavier than the blades. 

The most important dynamic phenomenon requiring verification by flight 





Vest is Made flutter. In order to study this phenomer. .1 under realistic 
conditions, it will be required that the frequency ratios and damping 
levels of the lower modes be approximately correct. Under these condi- 
tions, a scaling parameter which approximately measures the susceptibility 
to blade flutter is 

4/3 

( 1 ) 




where R = blade radius 


R g = distance to the sun 


T o = = the tension at the blade root 

M. = mass of one blade 


ft = spin rate (rad/sec) 
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It will be noted that the susceptibility to flutter is increased either 
by decreasing the spin rate, or by decreasing the distance to the sun. 
Thus, the severe conditions which exist near the sun can be simulated 
near the earth by decreasing the spin rate. Furthermore, since the 
build-up of unstable motions is very slow, flutter margins can be posi- 
tively and safely identified by first decreasing the spin rate until 
flutter occurs and then increasing the spin rate before the amplitude 
of oscillations becomes large. 

3.0 DESIGN CONSTRAINTS AND SELECTION OF PARAMETERS 

Figure 1 shows a comparison of parameters for the full-scale Halley 
vehicle and for a proposed flight test model. The parameters identified 
by asterisks (*) were selected on principles other than dynamic similitude 
design optimization, and may be considered to be design constraints. 
These parameters include the identity of the launch vehicle, the reflec- 
tive area, the film material, the closest approach to the sun, and the 
mission duration. The Ariane launch vehicle was chosen because of its 
availability and capacity. The selected values of reflective area, 
closest approach to the sun, and mission duration are based on the idea 
that the vehicle be capable of useful missions at the completion of the 
flight test program. The film material was selected to be the same as 
that for the full-scale vehicle. 

The blade chord (or more precisely, the deployment reel length) is 
a free parameter which was varied parametrically. The data shown in 
Table 1 correspond to a blade chord of 4.1 meters. Results for oti'.rr 
blade chord lengths are discussed in Section S. 


The remaining parameters in Table 1 were derived by scaling the 
full-scale Halley vehicle, using principles of dynamic similitude. 

Details of the scaling procedure are explained in the Appendix to this 
report. 

4.0 DESCRIPTION OF THE 4.1 METER FLIGHT TEST MODEL 

Figure 1 shows the 4.1 meter flight test model in its stowed con- 
figuration, and Figure 2 shows one blade in its deployed configuration. 
The vehicle has six blades which deploy in a single plane. The arrange- 
ment of the mechanical components, and also the deployment sequence, is 
identical to that of the full-scale vehicle, except that the central 
support column remains attached to the launch vehicle. Spaces are pro- 
vided between the deployment reels to permit bracing of the central sup- 
port column, if necessary to limit deflections during launch. 

The mass distribution of the 4.1 meter flight test vehicle is com- 
pared with that of the full-scale vehicle in Table 2. It will be noted 
that the mass fraction for structure is significantly smaller than that 
for the full-scale vehicle, in accordance with the square-cube law for 
the growth of structural weight with size. The mass fraction for elec- 
trical components is significantly larger, due mainly to the fixed size 
of many of the components in the control system. Detailed mass distribu- 
tions for the flight test model and for the full-scale vehicle are shown 
in Table 3. On balance, the mass per unit reflective area is slightly 
less for the flight test vehicle than for the full-scale vehicle. 

The total mass of the flight test vehicle (323.3 Kg) is small enough 
to permit large payloads when the Ariane is used as the booster. Figure 
3 shows the range of characteristic accelerations (acceleration due to 


1 
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normal incidence of radiation pressure at 1 A.U.) that are available for a 
range of payload mass within the lift capability of the Ariane. 

5.0 PARAMETRIC VARIATION OF BLADE CHORD 

The blade chord (deployment reel length) was varied from a minimum 
of 2.9 meters to a maximum of 4.4 meters. The upper value is the largest 
that will fit within the Ariane's payload envelope. The lower value cor- 
responds to a point where total mass is rising rapidly (see Figure 4). 
Detailed mass distributions for deployment reel lengths of 2.9, 3.4, 4.1, 
and 4.4 meters are tabulated in Table 3. Although the 4.4 meter chord 
produces the least mass, the 4.1 meter chord produces only a slightly 
larger mass and has significantly larger clearances with respect to the 
Ariane shroud. 

The large structural mass of the model with a 2.9 meter chord is 
mainly due to the higher blade root tension required to compensate for 
the increased blade radius, while satisfying the flutter criterion, 
Equation 1. 

i 

6.0 CONCLUDING REMARKS 



The Heliogyro model described in this report is capable of meeting 
the objectives of a useful flight test for the design of a full-scale 
vehicle of the size proposed for the Halley mission, and also for per- 
forming useful planetary missions, both within the constraints imposed 
by the Ariane launch vehicle. Smaller models, which might produce useful 
flight test data, but which would not utilize the full capacity of the 
Ariane, have not been examined in detail. It is estimated, from an 
earlier unreported study, that the minimum size vehicle which would be 
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useful as a flight test model has approximately 5000 m of refle rive 
area and a chord of one meter. This vehicle could easily be boosted to 
the minimum required circular orbit (1200 Km) by the Scout launch vehicle. 


Table 1 



o 


Comparison of Parameters for the Full-Scale 
Halley Vehicle and the 4.1 Meter Flight Test Model 


Parameter 

Full-Scale 
Halley Vehicle 

4.1 Meter 
Flight Test Model 

Launch Vehicle* 

Space Shuttle+I .U.S. 

Ariane (3-stage) 

Reflective Area* 

624,800 in 2 

60,000 m 2 

Film Material* 

.08 mil Kapton 

.08 mil Kapton 

Number of Blades 

12 

6 

Deployment Reel Length 

8.0 m 

4.1m 

Blade Radius 

7500 m 

2838 m 

Blade Aspect Ratio 

937.5 

692.2 

Closest Approach to Sun* 

.25 AU 

.30 AU 

Mission Duration* 

4 yrs 

5 yrs 

Rotational Speed 

.027 rad/sec 

.0514 rad/sec 

Blade Root Tension 

750 N. 

161 N. 

Relative Maneuver Time 

1.0 

.823 

Flap Hinge Offset/ 
Blade Radius 

.005 

.00587 

Maximum Flapping Moment 

1222 Nm 

56 Nm 

Design Torque for Pitch 
Motor 

5 Nm 

0.4 Nm 

Mass of Sail Module 

3837 Kg 

323 Kg 

Sail Module Mass /Are a 

2 

6.14 gm/m 

2 

5.39 gm/m 


* 

Design Constraints 


CP 
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Comparative Mass Distributions for the Full-Scale 
Halley Vehicle and the 4.1 Meter Flight Test Model 


Item 

Full- 

Halley 

-Scale 

Vehicle 

4.1m Flight 
Test Model 


& 

% 

k 

0 . 

'0 

Film + Coating ♦ Seans 

2123 

55.3 

208.4 

64.5 

Structure 

1436 

37.4 

77.0 

23.8 

Hardware and Mechanisms 

20Z 

5.4 

19.8 

6.1 

Electrical Components 

73 

1.9 

18.1 

5.6 

Total 

3837 

100.0 

323.3 

100.0 


Table 3 



Detailed Mass Distributions 


" " 1 - " " 

Deployment Reel Length 


.ieht Test Model: 


Full-Scale 
8.0 m 

. K F. 


2.9 m 

3.4 m 

| 4. 1 m 
iss (Kg) 

4.4m 


— ■ “ Me 


"" ' 

BLADES 






Bare Film (.08 mil Kapton) 

173.22 

173.22 

173.22 

173.22 

1764 

Coatings 

21.48 

21.48 

21.48 

21.48 

219 

Seams 

. 13.74 

13.74 

13.74 

13.74 

140 

Edge Tendons 

54.82 

40.56 

28.93 

25.56 

514 

Battens 

8.87 

8.87 

8.87 

8.87 

90 

Tip Weights 

3.00 

3.60 

4.44 

4.80 

20 

Hardware 

6.93 

6.93 

6.93 

6.93 

71 

Total 

282.06 

268.40 

257.61 

254.60 

2818 

RETENTION SYSTEM 






Flap Hinge Stays 

10. 10 

6.83 

2.57 

1.94 

61 

Deployment Reels 

14.00 

11.00 

10.70 

10.55 

145 

Yokes 

11.10 

6.67 

8.15 

9.47 

217 

Booms 

8.06 

4.43 

5.33 

5.93 

102 

Dampers 

.23 

.34 

.27 

.22 

11 

Pitch Axis Structure 

9.17 

8.12 

7.06 

6.70 

110 

Total 

52.66 

37.39 

34.08 

34.81 

646 

OTHER COMPONENTS 






Center Trv.ss 

5.39 

5.39 

5.39 

5.39 

197 

Pitch Motors 

6.92 

6.91 

6.84 

6.81 

48 

Pitch Bearings 

1.58 

1.40 

1.22 

1.16 

19 

Deployment Mechanisms 

6.93 

6.93 

6.93 

6.93 

84 

Control System 

10.00 

10.00 

10.00 

10.00 

15 

Electrical Wiring 

1.25 

1.25 

1.25 

1.25 

10 

Total 

32.07 

31.88 

31.63 

31.54 

373 

Total For Sail Module 

366.79 

337.67 

323.32 

320.95 

3837 



l 


i 
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APPENDIX 


DESIGN METHOD 

As stated in Section 3 of the main text, certain parameters were 
taken as given design constraints. These parameters are indicated by 
asterisks (*) in Table 1. The remaining parameters in Table 1 and also 
the detailed mass distribution in Table 3 were obtained by application 
of design principles. The principles and assumptions which were used 
to obtain each of the parameter and mass values are explained below. 

Where numbers are quoted, they refer to the 4.1 meter flight test model. 

1. Number of Blades 

The only number of blades given serious consideration was six. 

This is the smallest number of blades for which a symmetric (isotropic) 
rotor can be built which has damping in all modes. It is also attractive 
from the viewpoint of space utilization in the stowed configuration (see 
Figure 1). 

2. Sheet Width and Blade Radius 

The width of the reflective sheet was taken to be 0.4 m less than 
the deployment reel length in order ♦o allow room for edge tendons and 
gaps between the tendons and the sheet. The value 0.4 m is arbitrary 
and will allow for some flexibility in the design of the blade panels. 

As explained in the main text, the deployment reel length was treated as 
a free parameter. 

The blade radius was calculated as follows: 

R - (1) 


A-l 


2 

where A = reflective area (60,000 m ) 

= number of Hades (6) 

C = sheet width (3.7 m) 

The factor 1.05 is the allowance for cut-outs. 

3. Deployment Reel Diameter 

The considerations affecting selection of the inner and outer reel 
diameters are: 

a. Available space in the stowed configuration (see Figure 1) 

b. Adequate volume to stow the blade 

c. Adequate strength and stiffness 

d. Minimum weight 

The design formula for adequate volume is 

h < D o - D i> > 1 < 2 > . 

where D = outer diameter (.5 is) 
o 

D^ = inner diameter (.35 m) 

R » blade radius (2838 m) 

t = blade thickness allowance in the stowed configuration 
(.0283 mm) 

The value for t is the value used in the full-scale design. In the 
case of the 2.9 m design, it was necessary to reduce D^ to 0.3 m in order 
to satisfy Equation 2. 

4. Blade Root Tension 

The only consideration used in selecting blade root tension was 
that the susceptibility to flutter be the same as for the full-scale 
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vehic’i. An approximate design formula for equal flutter susceptibility is 




where T = the blade root tension 
o 

R = the blade radius 

r g = minimum distance from the sun 

The value of T q for the 4.1 meter flight test model (161 Newtons) 
is obtained from Equation 3 and the parameters listed in Table 1. 

5. Edge Tendon Hass 

The edge tendons must carry the blade tension while withstanding 
persistent attack by micrometeoroids. A i<*tional design procedure to 
meet the micrometeoroid threat was not available during preliminary 
design of the full-scale vehicle. However, in order to scale the value 
of edge tendon mass to the flight test model, it will be assumed that 
the cross-sectional area can be separated into a part required to carry 
ultimate tensile loads and a part that may be removed by micrometeoroids. 
It will be further assumed that each tendon consists of two .001 mil 
graphite polyimide tapes with frequent load transfer points. Under these 
assumptions, an approximate design formula for the mass of the edge 
tendons is 



where p = density of material (1527.5 kg/m 3 ) 

8 2 

* ultimate tensile strength of tapes (8.96 x 10 N/m ) 
= mission duration (5 years) 
and n^, R and T q have been defined previously. 
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When this formula is applied to the full-scale vehicle, it gives a 
value of M equal to 367 Kg rather than the value of 514 Kg given in 
Table 3. The value in Table 3 is considered to be overly conservative. 

6. Batten Mass 

It is assumed that batten weight is 5.1% of bare film weight, inde- 
pendent of vehicle size. 

7. Tip Mass 

Tip mass is assumed equal to 0.2 Kg per meter of sheet width, which 
is consistent with the full-scale design. 


8. Blade Hardware 


The mass of blade hardware is assumed equal to 4% of bare film weight 
independent of vehicle size. 

9. Mass of Film, Coatings and Seams 

The masses of these items are scaled directly from the full-scale 
design in proportion to reflective area. 


10. Spin Rate 

An approximate design formula for spin rate is 


ft 



(- 


.0514 rad/sec) 


where T q = blade root tension (161 N) 

= mass of one blade (42.935 Kg) 
R = blade radius (2838 m) 


(5) 


Equation 5 assumes that the blade mass is uniformly distributed 
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along the blade. 


— — 


1 1 . Relati v e Maneuver Time 

Relative maneuver time is defined as the time required to perform a 
maneuver with given blade pitch settings, divided by the orbital period 
around the sun. It can be shown that 

(T ‘C/nr,. 1 ' 2 *!) = constant (6) 


where T is the relative maneuver time, is the spin rate, r is the 
rm s 

distance to the sun, is the mass, of a blade, and C is the blade chord. 
12. Flap Hinge Offset 

Flap hinge offset is calculated to give the same value of blade 
flapping as the full-scale vehicle for a given error in blade pitch angle. 
It can be shown that 


r 2 RC 
f* s 


( 7 ) 


where e is the flap hinge offset and the other quantities have been pre- 
viously defined. For the 3.4 m and 2.9 m blades, the value of e/R given 
by Equation 7 eads to excessively large retention system weight and has 
been increased as shown in the following table: 


Model 

2.9 m 
3.4 m 
4.1 m 


e/r 
.006 ) 

QQ 5 f increased over values given by Eq. 7 
.00587 


4.4 m .00651 

Full-scale vehicle .005 
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13. Damper Stiffness 

The damper is designed to give the same damping coefficient for 
blade vibration modes as in the full-scale design. The design formula 
for damper stiffness is 


K 


d 


g l 


3T R 2 

-j- \i) (= 7086 Nm/rad ) 


( 8 ) 


where g^ = damping coefficient of damper (0.3) 

g^ = damping of first flapping mode (.001) 

14. Flap Hinge Stays 

The length of the flap hinge stays is computed by the following 
formula which maximizes flap hinge brace stiffness for a given mass. 





where K^ 
*b 


* damper stiffness 

= boom stiffness (assumed equal to 6.39 Kj) 
» yoke stiffness (assumed equal to 4.0 K^) 


The mass of the stays is computed by the formula 

16pT S’n, 

m x (= 2.57 Kg) 


Eh 


11 .2 

where E = modulus of elasticity (1.2065 x 10 N/m ) 
p = density of material (1522.5 Kg/m 8 ) 
h ■ vertical separation between stays (3.4 m) 


( 9 ) 


( 10 ) 
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In order to minimize weight, the vertical separation w r~, increased 
for the shorter deployment reels as follows: 


Deployment reel length h 

2.9 m 4.0 m 
3.4m 3.6m 
4.1m 3.4m 
4.4m 3.3m 


15. Mass of Deployment Reels 

The mass fraction of the deployment reels (relative to film mass) 
was assumed to be proportional to the square root of the deployment reel 
length. In addition, the masses of the shorter deployment reels were 
increased arbitrarily to account for increased vertical stay separation 
(see Section 14 above) and for decreased reel diameter (in the case of 
the 2.9 meter reel). 

16. Booms and Yokes 

The booms and yokes are designed by stiffness. The assumed ratio 
of boom stiffness and yoke stiffness to damper stiffness is given just 
below Equation 9 above. The cross-sectional area of effective material 
in the boom and in the yoke to provide the required stiffnesses is given 
by the design formula 





+ v 

Ew 2 



where i = length of yoke (4.1 m) 

= length of boom (2.68 m) 
w * width of cross-section (0.3 m) 

E * modulus of elasticity (1.2065 x 10 11 N/m 2 ) 


(ID 
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The masses of the booms and yokes are computed from 


^ = Sr^pJL^ (= 5.33 Kg) (12 

M y » 5n b p£ y A y (- 8.15 Kg) (13 

The factor of 5 is a structural inefficiency factor derived from the 
full-scale design. 


17. Maximum Flapping Moment 


The maximum flapwise bending moment is assumed to occur in the 
maximum precessional maneuver with one blade feathered due to a non- 
functioning pitch motor. It can be shown that the moment has the follow- 
ing proportion to design parameters: 


M ~ ( £_) C 

max V R / 


Table 1 shows that the value for the full-scale design is 1222 Nm, 
which gives 56 Nra for the 4.1 meter model. 


18. Damper Mass 


The design formula for damper mass is 


2n. pM 2 
b max 

Y 2 GK, 


(- .269 Kg) 


where n^ * number of blades (6) 

p * density of damper material (1384 Kg/m"*) 

Y « maximum shear strain (0.2) 

5 2 

G = shear modulus (6.894 x 10 N/m ) 

K^ » damper stiffness (7086 Nm/rad) 

M * maximum flapwise bending moment (56 Nm) 

insLX 



19. Pitch Motor Design 


The maximum pitch motor torque consists of a part due to the maximum 
static blade restoring moment and a part due to friction. The design 
formula is 


M e * * 675 x 750 x 


Gro) 


/7500\ 
X V R ) 


2C r r, M 
f b max 




.288 Nm) (16) 


where = coefficient of friction (.01) 

r^ = inner radius of bearings (.025 m) 

= distance between bearings (.15 m) 

M = maximum flapwise bending moment (56 Nm) 
max 

The factor .675 is the maximum blade restoring moment for the full-scale 
design. The design torque was increased to 0.4 Nm for conservatism. 

Tile maximum shaft power is assumed proportional to (Nq where ft 
is the spin rate of the vehicle. With this assumption, the power rating of 
the pitch motor for the flight-test model is 15.2% as large as that for 
the full-scale design. 


The masses of the pitch motor, of the pitch bearings, and of the 
pitch axis structure were scaled from full-scale design values by the 
following assumptions: 

2/3 

M (pitch motors) — (Power) (= 6.84 Kg total) 

M (pitch bearings) — (M ) 7 ^ 3 (* 1.22 Kg total) 

nicix 

2/3 

M (pitch axis structure) — (M ) ' (= 7.0b Kg total) 

IllttA 

20. Center Truss 


The mass of the center truss (5.39 Kg) was scaled from the mass of 
the hub end assemblies for the full-scale vehicle (90 Kg) by assuming 
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that the mass is proportional to the film mass multiplied by the two- 
thirds power of the distance from the centerline to the axis of the 
deployment reel in the stowed configuration. 

21. Deployment Mechanisms 

The mass of the deployment mechanisms (6.93 Kg) was assumed equal to 
4% of the bare film mass. 

22. Control System and Electrical Wiring 

The mass of the control system (10 Kg) was assumed equal to two- 
thirds of the mass of the full-scale control system. 

The mass of the electrical wiring (1.25 Kg) was assumed equal to 
one-eighth of the mass of electrical wiring for the full-scale vehicle. 
This was done because there are half as many blades, the distances are 
one-half as large, and the power rating is much smaller. 
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1.0 INTRODUCTION 


A rotating elastic body requires a dynamic analysis to determine the 
presence of dynamic instabilities. As part of the investigation of the 
Heliogyro blade design studies, a stability analysis was conducted using 
MSC/NASTRAN. This analysis required the development of solution methods 
for NASTRAN that included the following items: 

1. Centrifugal forces 

2. Tennis-racket moments 

3. Influence coefficients for blade flutter 

4. Structural damping of elastomeric damper 

5. Coriolis forces 

6. Hub boundary conditions for reac.tionle.ss , collective and cyclic 
modes analysis 

7. Blade pitch effects 

8. Vertical and inplane solar pressure loads 

9. Blade fabrication distortion effects 

The methods for implementation of the above items and the usage of NASTRAN 
for the stability analysis are the subjects for this report, including the 
unusual steps associated with obtaining a NASTRAN solution. 

Usage of NASTRAN to simulate the equations of motion for the Helio- 
gyro blades requires the addition of special terms which account for items 
1 through 5 above. Section 2.3 of Reference 1 presents a general descrip- 
tion of the equations of mot: on for a rotating blade. That description 
generally applies to the methods of solution developed for NASfRAN. Items 
1,2,3 and 4 were included with the stiffness matrix. Item 5 was included 
in the damping matrix. Item 6 was incorporated with the mass matrix. 
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Items 7, 8 and 9 were included through redefinition of the geometry of the 
basic NASTRAN model. Extensive alters to the NASTRAN rigid formats were 
necessary to generate the sequence of NASTRAN solutions including the 
effects due to items 1 through 5. 

2.0 NASTRAN RUN SEQUENCE DESCRIPTION 

This section describes the procedures that were used in the NASTRAN 
dynamic analysis of th'' Heliogyro Solar Sail. The intent of this section 
is to provide the reader (who is a NASTRAN user) with sufficient details 
to enable him to perform a similar type of analysis. The NASTRAN user is 
assumed to have some knowledge of the following NASTRAN capabilities and 
solutions: 

1. DMAP and Rigid Format Alters 

2. Differential Stiffness Solutions 

3. Modes Analysis (Direct and Modal) 

a. Real (SQL 3 or SQL 25) 

b. Complex (SQL 7 and SQL 10) 

4. Data Base Operations 

All of the above items are described in Reference 2. In addition to the 
above, it is important to understand the contents of Chapters 1 through 3 
in Reference 1. A description of the Heliogyro NASTRAN model and results 
of the analysis are described in Reference 3. 

2. 1 Straight Blades 

The designation of straight blades means that the blade is a straight 
line and the blade axis is coincident with the radial line emanating from 
the center of rotation. Only items which would cause no inplane or vertical 
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deflections ir.ay be analyzed with the following sequence of runs. The items 
included are centrifugal loads and/or imposed blade pitch. The sequence 
of NASTRAN runs which are required to obtain straight blade results is 
described by the following steps: 

1. The Differential Stiffness Solution (S0L 4,1) uses RF 4$33A 

which puts in the centrifugal stiffness terms and the tennis 

racket effect. The output to be used for the next run is the 

differential stiffness matrix [K ] which has been stored on a 

gg 

Data Base. When preparing the NASTRAN input, the following 
details should be considered: 

a. Set parameters M0DEL and THETAZ for storage of [K^ ] 

g g 

and the geometry data on the Data Base for doing structure 
plots. 

b. Select the proper set of DMIG data for the 0M2 (ft 2 data). 

The 2,2 term of 0M2 should be set to zero, corresponding 
to a prestretched blade. 

c. Select proper TENPAD data for correct tennis racket moments. 
See Appendix I for calculation of TENPAD. 

d. Eliminate unnecessary data such as ASET, EIGR, DMIG's 
for complex constraints. 

e. Select the proper set of SPC's. 

f. If this is the f'rst run and the Data Base has not been 
created, the following alter is necessary: 

ALTER 2 

DBMGR //0 $ 

g. Note the core allocations set with PARAM statements with 
parameter D0WN. 

h. Note whether EP0INT's or SP0INT's are included for later 
modes runs for the complex constraints. 
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i. When design changes for * v e blade are being made which 
include chord reductions, be careful with the RBE2's, as 
they don't have differential stiffness capability. 


2. React i onless Modes Solution (S0L 25) uses the alter RCTM0DE 
to fetch the differential stiffness matrix [K^ 1 stored on the 
Data Base during the S0L 4 run. The alter includes the centri- 
fugal force terms and tennis racket moments. The alter also 
stores the mode shapes for plotting at a later time. When 
preparing the NASTRAN input, watch for the following: 


a. Set parameters M0DEE and THETAZ for fetching the right 
[K^ ] from the Data Base. Set the parameter S0LID to 
identify the modes to be stored on the Data Base. 

b. Select the ASET for the Reactionless Modes. 

c. Remove the extraneous data used by the differential 
stiffness such as the SPC and RF0RCE data. 

d. Select the proper ASET, SPC, 0M2, C0UPLX and EIGR data 
for the type of modes. 


3. Collective Modes Solution (S0L 25) uses the alter RCTM0DE to 

fetch [K^ ] stored on the Data Base during the differential 
SS 

stiffness solution. This solution is nearly the same as the 
reactionless modes solution in item 2, except for the following 


a. Select the proper SPC's 

b. Set parameter S0LID for mode storage on the Data Base. 

4 . Cyclic (Regressive and Progressive Sequence) Modes (S0L 7,1 
or S0L 10,1) use alters NEWCYC or NEWFLUT, respectively, to 
put in the [K^ ] from S0L 4,0, the centrifugal force terms, 
the tennis racket moments and Coriolis force terms. When the 
"direct" approach is used in NASTRAN, structural damping is 
automatically included, however, the "modal" approach requires 
an alter to include the structural damping. The modal method 
was not the original approach used to solve for cyclic modes. 
The direct method was first used and only the inplane modes 
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1 


were obtained because of the core size limitation for the ASET 
points used in the analysis. Whenever preparing a cyclic modes 
solution, consider the following: 

a. Set parameters M0DEL and THETAZ for fetching the correct 

] for the run. 
gg 

b. Remove differential stiffness related data, SPC's and 
RF0RCE. 

c. Put in structural damping on CELASl's for S0L 7,1 and on 
CELAS1 and MAT1 cards for S0L 10,1. 

d. Select proper ASET, SPC, 0M2 and C0UPLX data. S0L 7,1 
uses SP0INTS for C0UPLX while S0L 10,1 uses EP0INTS. 

e. If doing a two-bladed model, special temporary storage is 
required due to the problem size. 

2 . 2 Blades with Inpl a ne and Vertical Deflections 

Whenever inplane or vertical blade deflections are included so that 
the blade axis does not coincide with the y-axis, the blade geometry 
changes are necessary to calculate the modes. Inplane and vertical 
blade deflections are caused by solar pressure loads and differential 
strain (manufacturing deformations) effects. The combined effects of 
solar pressure loads, differential strain, centrifugal loads and blade 
pitch may be analyzed in any combination with the following sequence of 
NASTRAN runs. 

1 . 


2. 


Differential Stiffness Solution (S0L 4,0) has to be performed 

to obtain [K^ ] for the straight blade (see Section 2.1). If this 
d *• 

[K ] was already obtained from a previous run, then this step 
g g 

may be omitted. 


Static Solution (S0L 24,1) uses the [ K gg] generated in Step 1 
and the STATS0L alter. The alter includes the addition of the 

straight blade [K^ 1 to the conventional stiffness matrix, the 

gg 
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centrifugal force effects and the tennis-racket moments. The 
solar pressure loads (see Appendix L) are put in to the bulk 
data by the conventional NASTRAN static solution through 
F0RCE data cards. Manufacturing deformations are introduced 
through the DEF0RM bulk data cards. The values put on the 
DEF0RM cards are 1/20000 times the rod element length on the 
trailing edge stiffener and +1/20000 for the leading edge. 

Solar pressure loads which simulate untrimmed collective pitch 
are shown in Appendix L. When setting up the NASTRAN deck, 
note the following details. 

a. Set the parameters M0DEL and THETAZ to fetch the proper 

[K^ ] from the Data Base, 
gg 

b. Select the proper loads, SPC's and 0M2 data (term 2,2 is 
-ft 2 now) . 

c. Remove the data associated with the modes solution for 
ASET, EIGR, SPC and C0UPLX. 

d. Remember to include the differential strains (DEF0RM) 
or inplane solar pressure loads to simulate untrimmed 
collective pitch. 

3. The Differential Stiffness Solution (S0L 4,0) has to be re- 
peated with the selected effects (centrifugal force, solar 
pressure loads, differential strain, etc.) included. In 
addition, the blade geometry has to be modified to include the 
deflections of the grid points obtained from the static solution 
(Item 2). Going from the static solution to this second dif- 
ferential stiffness solution has the effect of neglecting any 
internal element loads introduced by the static solution. To 
determine how well the >olar pressure loads are balanced by 
the centrifugal loads, vertical SPC's are introduced along the 
blade in this differential stiffness solution so that forces 
of single point constraints may be obtained to measure any 
unbalanced forces. The truss tip usually had significant SPC 
forces. Since internal loads were lost, primarily in the 
truss, the DEF0RM capability was used to introduce the internal 
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loads in the rod elements of the truss. The loads resulting 
from the DEF0RM cards were spread to the boom, yoke and post 
elements so that iterations were required to obtain a balance 
between the solar pressure and centrifugal loads. Details to 
check while assembling the NASTRAN deck are the same as the 
straight blade differential stiffness solution. 

4. The Modes Solution (S0L 10,1) uses the NEWFLUT alter and the 
[K^ ] matrix from the previous solution (Item 3). Since the 
blade has vertical and inplane displacements, the Coriolis 
forces must be included in reactionless, collective and cyclic- 
type modes solutions. The NEWFLUT alter includes the Coriolis 
forces, the differential stiffness matrix, centrifugal force 
terms, tennis racket moments and structural damping. Details 
of the NASTRAN deck to check are: 

a. Set the parameters M0DEL and THETAZ to get the proper 

[K^ ] matrix, 
gg 

b. Select the proper ASET and SPC data for type of modes, 
reactionless, collective or cyclic. 

c. Select the proper 0M2 data (term 2,2 is ft 2 ). 

d. Remove differential stiffness data, RF0RCE. 

e. Include structural damping on CELAS1 and MAT1 data cards. 

f. If cyclic modes, put in C0UPLX and EP0INT data. 

g. If two-bladed model, special temporary storage is required. 

2.3 Flutter Solution 


To obtain the aerodynamic force coefficients, calculate off-diagonal 
stiffness terms as described in Appendix M. The flutter solution is 
precisely the same as the modes solution (Section 2.2, Item 4) when 
vertical and inplane deflections are present, except for the addition of 
the AER0 term to be added to the stiffness matrix. 
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APPENDIX A 


S0L 4,0 Sample Deck for the Two-Bladed Model 


(RUN DBDS2B , 228 , DEAN , 25 , 300 /200 . ISO 

(a YM PRINTS,, HNS 
[ASG,UPR EBuprT(*i) 

(USE Z#£.HUPRT (♦( ) .. .. 

(BRkPT prints/* 

(QUAL ISDMSC 

(ASG, A . aNASTRAN, _ 

(HOG UPDATE RF9S33 

(ELT *RFALTEPS.RFaS33, TPFS.RFRS33 

(ASG, AX....DEANDB01* - 

(USE DU01,DEaNDB01 

(ASG, AX HIKEDATA, 

... (HOG . EDIT .0 A T A 

(ELT HlKEDATA.HigBLDl,TPFS,BLADEl 
(ELT HlKEDATA.MigBLD2 # TPF J.BLADE2 

. ,(E|,T_.HIKEPA.TA^M.15BLIQ*TPFS.BLD1T0 

(ELT HIKEDATA,H15B2T0,TPFS,BLD2T0 
(ELT HIKEDATA, HUB AQ07 , TPF S . HUBO A 

•1«,39 

•as, as 

*53,57 

_ (FREE_JIKEDATA, 

(ADD HIKEDATA. ASGCRDS 

NASTRAN TITUEOPT a 0 , DAYLIHIT a 99, wICORE a 90000 

ID HELIOGYRU SOLAR S A ILER_0 IFFEPENT I AL_ST I FFNESS RUN 

‘ TIME 6 
DIAG 8 
DIAC 22 
SOL 9,0 
ECHOOFF 

(ADD TPFS,RFas33 

S " PUT IN RF 'ALTER US33 

ECHOON 

CEND _ _ _ _ _ 

TITLE * HELIOGYRO solar SAILER 'two BLADED model 15# theta'z * 0 
SUBTI a OIF STF WITH CENT FORCE, Ra7500H# Ca?,63M,OMGa,027RAD/SEC 
ECHO a SORT __ 

SPC a 1 

ULQAD a ALL 
DISP a ALL 
SPCF a ALL 
elfowce a all 
GPFoRCE a ALL _ r 
SUBCASE’ 10' 

LABEL a STATIC SOLUTION, CUTOUTS AN 0 REDUCED CHORD 
LOAO a i 
SUBCASE 20 - 

LABEL « DIFFERENTIAL STIFFNESS SOLUTION 
DSCOEF a DEFAULT 


ORIGINAL PAGE 

of poor quality; 



BEGIN BULK 

HDD TPFS.BLADEl 

l ADD TPFS.BIA0E2 

[ADD * TPFS.BLDITO 

l ADD TPFS.BLD2T0 

[ADD TPFJ,MUBDA 

SPCl , 1,123*56, ROOOOO, 900001 
SPCl , 1 , 123*56, l 9997,29997 
P AR A M , MODEL ,15 
PARAN,THETAZ,0 

PARAM,GRDPNT,1 

enddata 
[BRKPT PRINTS 

[FREE ,H Z 

[SYM EBttPRT , , MN5 
[SYM PRINTS,, MNS 
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RF4$33A Alter for S0L 4,0 





BEGINNING OF R. ALTER «s33A 

♦♦♦♦♦♦♦♦♦♦♦♦♦ 23 NOV 1975 ♦ ♦♦♦♦♦♦♦♦♦* 

GRIDROINT FORCE BALANCE ANO ELEMENT STRAIN ENERGY RFl/2 

0a/i5/75 . _ RIGID F o R M A T a / SERIES Mi . . 

RF1/2 

CASE CONTROL INPUT RFi/2 

1 .REQUEST. FOR GRIDPOlNT BALANCED GRIDPCINTS SPECIFIED RF i/3 

THROUGH SET nsi, 2,,..- RFl/2 

RF1/3 

_ . GPFORCEsN RFL/2 

RFi/2 

OR RF1/3 

"GPFOR'cEsALL RFl/2 

RFi/2 

2 .REQUEST. FqR ELEMENT. S.TR A IN .ENERG Y_C ALCUL AT IUN5*FUR - ELEMENTS RFi/2 

specified through set rfi/2 

RFl/2 

ESE.SM RFi/3 

RFl/2 

OR RFl/2 

ESE’ALI RFL/2 


* RFl/2 

* — RFl/2 

ALTER 2 

PARAH //DOHN /25000 S 

DBMGR // S 

ALTER 8 

DBSTORE SIL,EQEXlN # ECT l CSTH,BGPDT//y,Y. M ODEL/V,Y,THETAZ S 


ALTER 25,26 RFl/2 

lM.fECT,EPT,8GPDT,SIL.GPtT,CSTM/ESTi.GEI J GpECt|/V 1 .N < LU$ET/C.Na23Z RFL/2 

V,N, N0SIhP/1/V,n,N0GENL/V,N,geNEL S RFl/2 

SAVE NOSIHP, nogENl,GENEL RFl/2 

CHKPNT E S T , G E I , G P E C t . S RF 1 / 2 

PARAM //AJD/V,N,NOKGGX/1/0 $ RFl/2 

PARAM //ADD/V,N,NOMGG/V,N,SKPMGG/l S UNTIL EMG Fix RFl/2 

ALTER 29,29 

CHKPNT OGPST J 
ALTER 30,35 

EMG __EST,CSTM,hPT 4 DIT,GEOm2 I4I /kEL m .KoiCT*.mELM,MDIC.T^,/.V.N,NOKGGX/ RFl/2 

V , N , NOMGG/0 ///C , Y , COUPMA58 S 

SAVE NOKGGX,NOMGG s rfi/2 

chkpnt kel m ,kdict,mel*,mdict. S 

EMA GPECT,KOICT,kELM,BGPdT,SIL,CSTM/kgGX,GPST S rfi/2 

chkpnt KGGX,GP3T % 

COND L^L 1 , SKPnGG _S 

COND JM>mgG,NOHGG RFl/2 

EMA GPECT,mDICT,mELM,8GP0T,SIL,CSTM/hgG,/C,N,.1/V,Y,i»TmasS » 1,0 S RFl/2 

CHKPNT MGG S 

LABEL JmpmgG * " RFl/2 


ALTER 97 

5 MULT 1 PLY TJmfs mgG, AQO fO KGG 
HATMOD OH 2 , , , , , / 0 M?G,/ 5 /V,N,LUSET S 

MTRXIN, , maTPOOL EQEXIN,sIL,/TENPAO,, /S,N,LUSET/S, n.noth s 

diagonal KGG/IGG/SQ'JAPE/O J> s 

ADD TENPAD, IGG/C 0 S 2 Th/(- 1 , 6 . 0 . 0 ) S 

MPYAD OM 2 G,COS 2 TH,/QMCnS S 



X 

1 


1 


J 

i 




MPYAD QMCnS,MGG,K3G/KGGQ S 

EQUIV KGGQ # KGG/AL^AYS 5 

ALTER 102 

PARAM //DONN/O S 

CPFDR CASECC,UGV,KELM,KDtCT,ECT,r:GEXTN,GPECT,PG , OG , BUPDT . S IL , CSTM/ RK1/3 
r)NRGYl,OGPFBl/DSO/C,Y,TINY S 
PARAM //DOWN/23000 S _ 

OF P UNRGY1 ,0GPFB1 , f ,, // S RFi/3 

ALTER 103,107 

CONO ... P 2 , JUMPP |_07 s 

PLOT PLTPAR,GPSETS,EL3ETS#CASECC,BGPDT,EQE2IN,SIL,PUGV1 , ,GPECT,0ES1/ RF1/3 
PLOTx2/V,N,NSIL/V,N,LUSET/V,N. JUMPPLOT/V,N,Pl ."/V,N, PFILE RF1/3 

SAVE PFILE S 

PRTMSG PLOTX2 // S 

LABEL P2 S 

ALTER 108,109 

P AR AML CASECC//0Ti/i/7//V,N,TSET S 

PARAML CASECC//0TI/l/6//V,M,0EF3ET S 

EMG EST, CSTM, MpT, DIT, , uGV , GPTT , ED T/ kDELM , KDD IC T , , , ,/ l/0/0/ 

///V,N,T3ET/V,N,0EFSET S . 

CHKPNT KDELM,kDOICT S 

EMA GPECT,KOOICT^KOELM,_8GPOT,SIL*CSTm/KOGG,/-1_S 

CHKPNT kDGG S 

06ST0RE K0GG//V, Y , MODEL/ V , Y , THETAZ/ J 
OBMGR // S 
'EXIT S 

ALTER 160, 1 b^> 

COND P3,JlJMPPLOT S 

PLOT Pi. TP AR , GPSE TS , ELSE TS # C ASECC , BGPOT » EGE? IN, S IL , PUBGvi , , GPECT , OE5B 1/ 


PLOTX3 /V,N, NS IL/V,N,LUSET/V,n. JUMPPLOT/V,N,PLTFLG/V,N, PFILE RFl/3 


PRTMSG PLOTX 3 // 1 
LABEL P3 3 

S ITERATED DIFFERENTIAL STIFFNESS JRF_a/ 

s io/ 1/73 rigio format a /series mi rf a/ 

s RF «/ 

S CASE CONTROL INPUT RF _fl/ 

S OSCOEFFICIENT CARO RF «/ 

s R 6 tt / 

1 BULK DATA INPljT RF a/ 

J OSFACT card RF «/ 

S PAHAM?S EPS IO, NT, BETA ARE OPTIONAL RF «/ 

I _ RF a/ 

1 IMPROVES accuracy of differential stie f ness MATRIX By ITERATION RF «/ 

••»---RF «/ 

ALTER 98,98 * _RFll^21 RF__«/ 

SDR1 USET', ,ULV#UO0v#VS,G0,GM,Ps,KFS,KSS,/UG'V,PG 1 ,OG/C,N,l/C,N, RF «/ 

OSO 3 RF a/ 

ALTER 99,100 RF a/ 

CHKPNT UGV,OG S ‘ " RF a/ 

SDR2 CASECC, CSTM, MPT, DIT, FQEXIN, STL, GPTT, EOT, BGPOT,PG,GG,UGV, EST,/ RF 0/ 

UPGl ,0GG1 ,0UGV1 ,f)ESl ,0EF1 ,PijGVt/C,N,CS0 S_ _ RF 0/ 

♦ITER 110'" " • " RF a/ 

PARAM //C,N, AD0/V,N, ShIFT/C, N, -1/C, N,0 1 RF «/ 

PARAM _ //C,N, AOD/V,N,COUNT/V,N, ALMayS 5-1/V,K,NEVER*1 S RF a/ 

PARAMR ■ "/7C,N,ADO/V,N,DsEPSI/C,N,0,0/C,N,e,0 S RF a/, 

PARAML YS//C,N, NULL/C, N, /C , N, /C , N, /$, N, M0YS S RF U/ 


S, YS,/P3S/C,N,0/C,N, 1/C,N, 1/C,N, 1 ! 

,ys,/pfs/c,n,o/c,n, i/c,n § t/c#n , 1 s 

T, PFS,PSS/PN/C-,.N,N/C,-N,F/C* W »S j 

PEX/mPTF 1 


JUMP OUTLPTOP 

LABEL OUTLPTOP S 

EQUIV PG,PG1/nOYS 3 

CHKPNT PCI . .. 

PARAM //C,N,KLOCK/V,N,TO 

ALTER 129,136 

ADD KAA,KDAA/KBU.-t — 

ADD KFS,KDF3/KBFS * 

ADD KSS,KDS3/KBS8 3 

CONO - PGOK r NOYS-J- 

MPYAD KBSS 

MPYAD KBFS 

UMERGE . USE! „ .... __ 

EQUIV PN,PGX/mPCF1 

CCMD NOMS,MPCFl 

UMERGE USET,PN r /PGX/C,N,G/C,N,N/C,N« 

LABEL NQMS s 

AOD PCX, PG/PGG/C,N(,1, 0,0,0) 

EQUIV PGG,PG1/ALRAYS 

LABEL PGOK 

AOD PG1,/PG0/ 3 

ALIER tat 

JUMP 
LABEL 

PARAM 

SSG2 
ALTER 

SDR 1 

ALTER 
ADD 

EMG ... ESI 


M_ S — 


INlP TOP 
IMLPTOP 
_//C,N,KLOeK/V,N^TI 

i t r r • r u u « » » r n r /« 


RF a/ 
RF a/ 
RF 4/ 
RF 4/ 
RF 4/ 
RF 4/ 
RF 4/ 
RF 4/ 

RF~4/ 

RF 4/ 
RF 4/ 

RF- 4/ 

RF 4/ 
RF 4/ 

RF— 4/ 

RF 4/ 
RF 4/ 

RF-4/ 

RF 4/ 
RF 4/ 
RF_4/ 


s 

USET,GM, Y3#KDFS,GO, ,PG1/»PBq.PBS,PBL S 
148,148 

....USEJT,.,UBLS#.,-YS # .G0*GM,PBS,KBFS*KBS3#VUBGV # ,.QBG*C r N,l/C,N, 
15T3, 156 

UBGV,UGV/DUGV/C,N, ,0,0,0 J S 

,CSTM*MPT,DIT,,DUGV,GPTTrEDT/DKOELM,DKDtHCT,,.,.,./l/a/0/ 

///V,N,Ts£T/V,N,oEFSET S 
CHKPNT DKDELM,DKDDICT S 

EMA . GPEtt*DKQDLCl^DKDELM > BGPDT^SIL*CstM/DKDGG^/ 

CHKPNT DkDGG S 


RF 4/ 

— RF_4V 

RF 4/ 
RF 4/ 
D31 -1- RF-4/ 


RF 4/ 


CHKPNT 
MPYAD 
DSCHK 

SAVE 
COND 
CUND 
EQUIV 
EQUIV 
EQUIV 
REPT 
TABPT_ 
LABEL 
ADD 

CHKPNT 

EQUIV 

CHKPNT 

EQUIV 

REPT 

TABPT 

label 
alter 

ALTER 
$ End 


N , 


DKDGG S 

DKDGG,UBG9,PGQ/PGI1/C,N,0/C,N, 1/C,N, 1/C,N, 1 S 
_PGL,PGIt ,UBGV//C#-Y,EPSIQsi r 5F-5/V^N,0SEPSI/C# Y,NTs. l/v^N 
TI/V,N,D0NE/V,N,SHlFT/V,N,CnuNT/C,Y,8ETAs4 S 
OSEPSI, DONE, SHIFT, COUNT S 

DONE, DONE 


DONE, DONE 

SHIFT, SHIFT S 
PG.PG1/NE7ER S 
PGlL.PGt/ALWAYS 3 
PG1,PGI1/NEVER S 
INLPTOP, l 000 S 
PGH,PGl,PG,»/Z _S 


SHIFT 3 

DKDGG,KOGG/KDGGl/C,N, (»1 ,0 ,o'.0) 3 
KDGG1 3 ... . . .... _ 

UBGV, IJGV/AL" A YS /KDGG 1 ,KDGG/ ALWAYS 
UGV,KDGG 3 

KDGG,KDGGl/NEVER/UGV,UBGV/NEvER s. 
OUTLPTOP, 1000 3 
KDGG1,KDGG,UGV,,// % 

_ . dune s 

170,171 

176,177 

OF RF ALTER. 4333 


tA0 * » 
w TO » quautx 



Q 
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S0L 24,1 Sample Deck for One-Bladed Model 



ISO 


(RUN 0bSTAT,228 , DEAN, 06, 300/200 
(SYM PRINTS, ,MNS 
(ASG, UPR EBlPWT(*t) 
tUSfc Z,EB1PRT(*1) .. .... 

(BRKPT PRINT S/Z 

• tQUAL ISOMSC 

(ASG, AX PEANOB0I, 

(USE ohoi,oeandboi 

(ASG, AX MIKEDATA, 

. (HOG _ EDIT. DATA 

(ELT HlKEDATA,MUBLOi,TPFS,0LAOEl 
(ELT MIKEDATA, M UB 1 TO , TFpJ , BLD 1 TO 

(ELT MIKEDATA, HUBA0D7*TPFS,HUBDA 

-1«,39 
*49,49 

•53,57 

•59,73 

(Elt mikf.oata,$olprs,tfps,solprs 

•1,31 __ 

(ELT MIKEDATA, STaTSOL, TPfS.STATSOL 
(FREE mikEDATA, 

(ADD. MIKEDATA, ASGCPOS 

NASTRAN TITlEOPT a 9 , DAyLIMIT * 99 
ID HELIOGYRQ SOLAR SAILER SOLAR PRESSURE LOADS 

TIME_2_ 

DIAG 8 
0 1 AG 22 

_SOL 24,1 

(ADD TPFS.STATSOL 
CEND 

TITLE s_HELIOGYRO solar_sailer one blabed MODEL i 0, THETAO * 0 
SUBTI s SOLAR PRF3SURE LOADS AND ONE DFG UNTRIMMED COLL PITCH 
VJ W LABEL = R = 7SOOM,La7,63M,OMGa,o27RAO/SEC, CUTOUTS ANO REDUCED CHORt 

ECHO = SORT 

SPC ■ 1 

OLOAD * ALL 

OISP * ALL 

SPCF = ALL 
ELFORCE a ALL 

GPFORCE = ALL 

LOAD s 2 
BEGIN BULK 

(ADD . TPF S, BLADE 1 

(ADD TPFj,BLD1T0 
[ADD TPFS.HUBDA 

(ADO _TPFS,30LPRS 

SPCl ,l , 123456, 900000. 19~997 
SPCl.l,, 999990 
PARAM, MODEL, 14 
PAWAM, THETAZ, 0 
PARAM, GROPNT, 1 
enodata 
(BHKPT PRINTS 
(FREE,R Z 
(SYM EB1PRT,,MNS 



s beginning of alter to get defections of rotating blade 

S *IT^ SOLAR pressure 

alter Q3 

DBFETCH /*DGG, , , ,/0,Y,MnpEL/V, Y. THETaz/ * 

ADD KGG, KDGG/KGGDS/ S 

O MATmOd 0M2,,,,,/0M2G,/5/V,N,LUSET S 

MTRXIN, , M ATPOOL,EQEXlN # slL, /TENPAD, ,/S, N,LUSET/S,N, NOTH S 

DIAGONAL KGG/IGG/SOUARF./ft t 0 S 

ADD TENPAD, lGG/COS2TH/(-l .0,0,0) S 

HP Y AD 0M2G,C0S2TH,/0MC0S % 

HP Y AD OHCOS, HGG# xGGDS/kGGQ S 
EQUIV KGGQ,KGG/ALHAYS $ 

alter ioa,io5 

* END CF ALTER 





ORIGINAL PAOB ® 

OP POOR quality 
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S0L 25,1 Sample Deck for Two-Bladed Model 


c* 




[RUN 0BHD2B.228 , DEAN, 25 , 300/20 0 , ISD 

IS YM PRINTS,, HNS 
{ASG.UPR EB5PRT(M) 

IUSE Z,EBSPRT(M) .... - . . . 

IBRKPT PRINTS/Z 
tOUAL ISOMSC 

(ASG , A *NASTRAN r . - . ..... .... . 

I ASG, A X DE ANOBO l # 

(USE DB01,DEANDB01 

(ASG , AX -MIKEDAT A r 

(HDG EDIT DATA 

CELT HIKEDATA,Mt5BLDl,TPFS # BLADEl 

IELT... MlKEDATA,Ml5BLD2^TPFS r BLADE2 : 

CELT *IKEDATa,m 15 biT0,TPFS,BLO1T0 
IELT MIKEDATA .MJ582T0, TPFS.BLD2T0 

CELT MlKEDATA r HUBA0D7,TPF$,HUBDA 

• 9,16 
•28,39 

-«9,<i9 

•53,57 

(ELT MIKEOATA,RCTHODE,TPFS,RCTMODE 

(FREE MIKEDATA, 

(DELETE, C M15CHD0ICT, 

(ASG, UP H15CMODICT,P///50 

(BR KPT. RUNG HS/mi5C^DD1CJ 

[DELETE, C H15CHDNPTP, 

(ASG, UP Ml5CMDNPTP,F///500 

(USE NPTP,H15CMDNPTP 

[ADD MIKEDATA, ASGCRDS 

NASTRAN TITLEOPT s 0 , OAyLIMIT * 99 , HICORE s 82000 

ID HEL10GYRO SOLAR S AILER-VIBR A.T.JON- MOOES RUN 

CHKPNT YES 
TIME 10 

DIAG e 

SOL 25,1 

(ADD TPFS.RCTMQDE 

TITLE s HELIOGYRO SOLAR SAILER TWO BLADED MuCEL 15, THETAZ * 0 
SUB T I s MODES WITH GENT FORCE, Rs7500m. C*7,63M, OHGs , 027R AD/SEC 

LABEL. .3 REOuCED .CHORD .AND .CUTOUTS^-COllECTLvE-MODES • 

ECHO s SORT 


SPC s 1 

METHOD * 2 . 

DISP * ALL 
BEGIN BULK 

(ADD ,_TPFS.BLADE1 

(ADD TPFS.BLADE2 
(ADD TPFS.BLD1T0 

(ADO TPFS.8LD2T0 . 

(ADD TPFS t HUBDA 

ASET1, 36, 900000, 9ft0001 

EIGH,2, M Glv,0 t . ,02,3a,3R,,U-9,*EIGR2. 

PARAM, SOLID, 3 

SPC1,1,12«5,90000(}»900001 

SPCl, 1, 12«56, 19997,29997 


ORIGINAL PAGE IB 
QF.J^&.QUALITfl 



PARAM, MODEL# 15 
PARAH, TKETAZ* 0 
PARAM, GRDPNT, 1 
ENDDATA 
IBRKPT PRINTS 
IF REE , R Z 
tSYH EB5PRT, , MNS 
tSYM PRINTS,, MNS 









APPENDIX F 


RCTM0DE Alter for S0L 25,1 




$ BEGINNING of alter 

s RIGID FORMAT 25 A^ TER TO GET RFACTIflNLESS MODES FOR 

J A STRAIGHT BLADE WITH A CENTRIFUGAL LOADING 

ALTER 92 .1 OIFF, 3 TFF , , CENT. FORCE Am C TENnIS RACKET HOHENT 
DBFETCH /KDGG, , , ,/9,¥,MQpEL/V # Y, THET a?/ S 
ADD' KGG,KDGG/KGGDS/ $ 

MATMOD OM2r i f 1 1 /OM2G* /5/V 4 N,LUSET.S. . 

MTRXIN, ,hatpqol,egexin,sil#/tenpad, ,/S,n,lusET/S.N # NUTH S 
DIAGONAL KGG/IGG/SQUARE/0,0 s 

ADD TENPAD, IGG/C032Th/(-1 .0«0.0)._S 

MPYAD 0M2G,C0S2TH,/0MCnS S 

MPYAD OHCOS,HGG,KGGDS/kGGQ J 

EQUIV _ KGGD, KGG/ALWAYS § _ 

ALTER 102, 103 S AvOID SOL ?"*f" ERROR 

ALTER 15U S SAVE BLADE MODE SHAPES FOR PLOTTING 

OB STORE PHIG,LAMa//v, Y, M0DEL/Vi.Y # SOLID. * 

DBMGR //7/V, Y, MODEL/V# y» sol id/v, Y# HqDEL/V# Yi SOLID/O/PHIG/UGV/ 

dbhgr // s 

S END OF ALTER 
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S0L 7,1 Sample Deck for One-Bladed Model 



- - • « ' 

TRUN,n C YCmoo ,228 ,DEaN,23,uoo/20o , ISO 

(SYM PRINTS,, HNS 

[ASC,UPR EU3PRT(*n 

tUSE Z , EB3PRT ( ♦ 1 ) .... 

(RRKPT PRINTS/Z 
IOUAL ISOMSC 

(ASG, A tNASTRAN, .. . .. . 

(ASG, AX MIKEDATA, 

{ASG , AX DEANOB01, 

tUSt. - D B 0 1 , OEANDBol 

(HOG EDIT DATA 

{ELT HlKEDATA,HUBLDl,TPrS t BLAOEl 

-261,28 2__ 

(ELT MlKEDATA,MlaBlTO,TPFS,8LDtTO 
(ELT MIKEOATA.HUBAOD7,TPFS,HUADA 

•10,2J 

-a9,«9 

(ELT MlKEDATA # NEwCYC,TPFs,NENCYC 

(FREE. HIKEOAT.A, 

(DELETE, C H1UB1CYDICT, 

(ASG, UP H1UB1CY0ICT,F///100 

(BRKPT PUNCHS/MUBIC.YOIC.T 

(DELETE, C MjaBlCYNPTP, 

(ASG, UP M1UB1CYNPTP,F///1000 

(USE_NPTP,HiaBlCYNP.TP 

(ADO MIKEDATA, ASGCRDS 

NASTRAN TITLEOPT a 0 , DAyLIMIT a 99 

±t . H±Ll 0GYRQ--S0LAR-3A I LER_ CYCLIC— MODES— RUN 

TIME 7 
CHKPNT YES 

. OIAG 7 

D I AG 8 
OIAG 13 

. OIAG 20 

'SOL 7,1 

(ADD TPFS.NEwCYC 

CENO 

TITLE = HELIOGYRO SOLAR SAILER ONE BLAOED MODEL U, THETAZ * 0 
SUBTI a CYC MOS MTH CENT, LOADS, Ra75flOM, Cs7,63m, UMGa , 027R AD/SEC 

LABEL *. C = 8M_ AT -l.<«5»,61R, CUTOUT- -AT.-4-,9«,95lR 

ECHO a SORT 
SPC a 1 

CMETHOD a i_. 

DISP a ALL 
BEGIN BULK 

(ADD ...TPFS.BLADEL 

(ADD TPFS.BLD1T0 
(ADO TPFS.HUHOA 

ASET, 999990, 0,900000, 12 

PARAM, MODEL, U 



PARAM,THETAZ,0 

CELASt, 19992, 19992, 19996,5,19997,5 
PELAS, 19992, 100, ,,3 
MATl,ll,120.65*9,,,3,, t , # 3 
ENDDATA 





(BRKPT PRINTS 
(FREE,R Z 

. (SYM EU3PRT , ,.MNS . 
(SYM PRINTS, ,MNS 




APPENDIX H 


0 

NEWCYC Alter for S0L 7,1 




RIGID FORMAT ALTER TO INCLUDE THE EFFECTS OF 


centrifugal forces 

CORIOLIS EFFECTS 
tennis racket effects 


a». COMPLEX constraints at. the. center of rotation .... 


rotating bodies. 


ALTER 33 S 

SETVAL _//V i N,NO0GG/O_S_OVER.RULE_SHA2 

ALTER U7 S 

MATMOD 0M2,siL,,,,/0H2G,/5/ S 

D 8 F E_T C H _./ 1 K D G G , , , , / 9 , Y « M 0 D E L /.V . Y , T.H E I A 2 . * 

ADD KGG,KDGC/KGGA 5 

EOUIV KGQA, KG G /ALWAYS S PAP 

HTRXIn, .^ATPOOL^ESEXIN.sIL^/TENW^jCcUPLX^S^^LUSET/SiN^OJH.J 

DIAGONAL KGG^IGG/SOUARE/O.O S 

AOD TENR-*£,lGG/COS2 T M/f-l. 0,0,0) S 

HPYAO CH2.G 4 C0S£TH/l/QHCQS..S 

MPYAD 0HC0S,hGG, /OHMGG S 

AOD OHMGG, /KCENT/ $ 

MATMOD OHt t SIL#,t t /OHlG,/5/_..S 

MPYAD QM i G # MGG , /BCOR 1 

ADD BCQR f BGG/BGGQ/ S 

EQUIV 0GGQ l BGGzALWA_1S_S 

ALTER 136 S 

UPARTN USETD.COUPLX/COUPLXD, ,, /P/D/O S 

maTWEDU_KCENT i USE_T,Gm,GQ/kaaQ/S,N,N0kcENT_S 

ADD KDD,KAAQ/KDOQO S 

EQUIV KDDQG,KDD/ALWAYS S 

ADD MOD j COUP L.XQ/MODQ S 

EQUIV hdoq',mdd/always i 

ALTER t«2 S GENERALIZED MASS CALCULATION 

MPYAD _ MOD, PHID^/MPHJ_S 

MATHfjD PHID,,,,,/PHIDCNJ,/10 S 

MRYAD PHIDCNJ,mPhI,/GMCX/1 S 

M-TMOD GMCX, ###t /GMCXF,/2////L fc -5 S 

MATPRN GMCXF // S 


ENO OF" RIGID FORMAT ALTER 
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S0L 10,1 Sample Deck for Two-Bladed Model 


O 


(RUN,N/T D6C M 2B,228 , DE AN , ? 3 , 2500 /<* 00 

[SVM PRINTS,, HNS 

tASG, UPR EBbPRT ( ♦ 1 J ,F///5000 

(USE Z»tB6PRT(tl) 

[BRkPT prints/z . 

(OUAL ISDMSC 
(ASG,A *NA5TRAN t 
(ASG,AX MIKEDATA, 

(ASG,AX DEANDH01, 

(USE DB01 ,DEANDB0t _ 

(HOG EDIT DATA 

(ELT MIKEDATA,M15BLD1,TPFS,BLADE1 

-281,282 

(ELT mIk£DaTa.mI'5HID2,TPfS i bLaDE2 . 

•278,278 

_ (ELT MlKEDATA t M15BlTO,TPES,BL01TO 

(ELT HIKED AT A,Mi5B2T0,TPFS,BLD2T0 
(ELT MIKEDATA,HUBA007,TPFS,HUBDA 

•9,1 6 

-28,31 ‘ 

-«9,«9 

-59,59 

EPOINT, 999990 
-73,73 

EPOINT, 999991 

(ELT mikedata, neVfl"ut",tpfs,cplxmd~ 

(FREE MIKEDATA. 

_( DELETE, C_H 15B1CPDICT, 

(ASG,UP Hl5BlCPDICT,F///300 
(BRKPT PUNCHS/M15B1CPDICT 
CSAVETAPE NPTP, 

H15-2BLD CH' NPTP 
(ADO MIKEDATA, SECONDS 

(ADD MIKEDATA, A_SGCRDS_ 

NASTRAN tlf'LE'OPT sO, DAyLIHIT s 99 , MJCOBE * 

ID meliogyro solar sailer complex modes pun 

TIME 21 
CHKPNT YES 
DIAG 19 
OIAG 7 
0 1 AG 8 
OIAG 13 

OIAG 20 

SOL 10# 1 

(ADD T?F JjCPLxMD 

CEND 

TITLE' s HE’. IQG YRO SOLAR SAILER TWO BLARED MODEL 
SUBTI s Rs7500M, Cs7,63M, OMGs’, 027R AD/SEC 
LABEL s REDUCED CHORD AND CUTOUTS, CYCLIC MODES 
ECHO s SORT 
SPC s l 
METHOD * 2 


127000 


MODEL 15, THETAZ * 0 
MODES 




I 


SET 100 



1, -7, 2, -7, 3, -7, a,-7, 5, -7, 6, .7,7, -7, 8,-7, 9, -7, 10, -7, 

It. -7*12. -7, 13. -7, m.. 7, 15,-7, 16. -.7, 17,-7, 18,-7, 19. ,7, 20,-7, 

2 1 .-7, 22.-7, 23,-7,29.. 7, 25,-7, 26.-7,27, -7, 28,-7, 29.-7, 30. -7, 

31 ,-7,32, -7, 33, -7. 3a, -7, 35, *7, 36, -7 ,38,-7,39,-7,40, -7 
SFWEQ = 100 
OISP = ALL 

BEGIN BULK . . . _ __ .. . 

I ADO TPFS . BLADE 1 
(ADD TPFS , BLAPE2 

IADD --TPFS , BLD 1 TO- 

l ADD TPFS.BLD2T0 
(AOD TPFS , HUBO A 

ASETL, 0,99999 0,99999 1 1 1 

ASE n , 12,900000,900001 

CELASi , 199 92, 1999 2, 1 9996 , «j , 1 9997 , 5 

CELASl, 29992, 19992, 29996, 5, 29997,5 

PELAS,19992,100,,',3 
M ATI, li, 120, 65*9,, ,3, ,,,,3 

t IGR , 2 ,MG I V , 0 , , , 02 »aO , 00,, t ,-9 ,*E I'3R2- 

PARAH,LHQDES,UO 

PARAM, model, 15 

PARAM.THE-TAZ.a 

PARAM, GRDPNT , 1 
ENDDATA 

[BHKPT. PRINTS 

IFREfc , R l 

tSYM EB6PRT, , MNS 

[S Y.H_PP INTI ,^HNS 



» 


APPENDIX J 


© 

NEWFLUT Alter for S0L 10,1 



0 


S BEGINNING GF alter for CQHPlEx MODES 

s R,F. 10 ALTER FOR COMPLEX MODES 

s modified for flutter mr is may 77 

% INPUT AERO FACTORS IN D m IG AERO* And SCALAH FACTOR IN 
s param acoef, complex 
alter a s 

PARAM //DQWN/53000..S . 

ALTER 2l,2i S EMG AND EMA For R,F, 10 

T A 1 , ,ECT,EPT,BGPDT,SlL,GPTT,CSTM/pST* , GE I , GPEC T , /V, N, LUSF. T/ 1 23/ 

Y tN, NOS IMP / l/V *_N # .NOGE.NL/-V*N» GFNEL ..i> 

ALTER 25,25 S 

CHKPNJ EST,GEI,r,PECT,OGPST 3 

SETVAL t/J, N, NOkGGX/1 $ | 

SETVAL //V,n # NC1mGG/1 S 

SETVAL //V,N,NQBGQ/l S 

SETVAL A/V,N,NO.K«GG/L.S 

EMG EST,CSTM,mpt,DIT,GE0M2, ,,/kElM,KDICT,MELM,mdICT,BELM, 

8DICT/S,N,NOKGGX/S,N # NOMGG/S # N,NOBGG/S,N,NOKOGG//C, y, 

coupmass s_ 

CHKPNT KELM # KDICT,MELM,MDICT,8ELM,BDICT S 

alter 26,2b s 

EMA GPEC T^KpLCT t^.ELM^BGPD.T^S LL*CSI*/KGGX*.GPSL_S 

ALTER 28,30 S 

COND errori,nomgg S 

EmA jPECTtMOlCt,MELM,BGPOT,3ILtCa.T M /MGGt/»l/V,y i wTMASSsUO _3 

ALTER 31 3 

PURGE BGG/NOBGG S 

COND lE_MA0 t _NOBGG_S 

EMA gpect^dict.b'elm.bgpot^il^stm/bgg^s^nobgg S 

LABEL LEMAB S 

CHKPNT 8GG 

PURGE Kar,G/N(]K/?GG S 

COND LEMAKa,NQK«GG S 

EMA GPEC T lK Q LC.T 1 K E LM^H G PQJ tSILi.CS tM /JL« G G X J S f jN t NO k aGG .3 

CHKPNT KaGG S 
LABEL LEM AKfi S 

ALTER «2 S CENT ,.AnD_C OR LOLI 3 -FORCES 

matmqd 0M2,,,,,/0M2G,/5/v,N,LUSET 3 
DBFETCH /KDGG,,,,/V,Y,MODEL/V,Y,TH£TAz S 

ADD _ kgg,kdgg/kgga_s 

CHKPNf KGGA S 

EQUI V KGGA , KGG/ ALWAYS 3 

CHKPNT KGG 3 

M T R X I N , , M A T P 00 L , E U E X I N , S tV , / T E N P A D , , / S , N , L 0 S E T / S’ , N , N 6 T H V 
DIAGONAL KGG/IGG/SOUARE/0,0 

ADD T E NPA 0 , 1 GG /.C 0 S 2 TM/fl _.0 ,.0 , 01_3 

MPYAD 0M2G,CnS2TH,/0MC0s $ 

MPYAD OMCOS,MGG, /QMMGC- 3 

ADD OMMGG, /KCENT/3 . . 

CHKPNT KCENT 3 

MATMUO OM 1 , , , , , /0 M 1 G , /5/V, N , LUSET 3 

MPYAD f)M 1 G, MGG, /BQQR 3 _ _ _ . 

CHKPNT BCOR 3 

S 1. 

s 2.. . _ SMP1...T.YPE. OPERATIONS. _ . 

alter Ru 3 print REAL mooes 

MATGPR GPL,USET,SIL,PHIA//H/A % 




; 


. 


ALTER 105 S 

MTRXlN, . M A TPUnL,EQDYN, ,/CDUPLX, AERO,/V. n,LUSETD/S,N,NOCOUPLX/S,N, 
NOAERO S P-SIZE MATRICES 

EQU1V Kacr, f KUPP/M0UE/KCENT#KCPP/NnUE/BGC # BPP/NCUE 
EOUIV BCOM,9CPP/NiOUE/Kr,G # KPP/MC)UE S 
CHKPNT KaPP,KCPP,BPP#BCPP # KPP S 

COND LNOPSET,NOUE * _ _ _ 

VEC USETD/VPGE/P/G/E S 
MERGE K4GG,,,,VPGP,/K4PP t 

MERGE KGG lllf VPGE,/KPP J 

MERGE KCENT, , , , VPCE,/KCPP S 
MERGE BGG, f ,,VPGE,/BPP S 

MERGE BCOR, , , , VPGE#/BCPP_ S _ 

LABELLNOPSET S 

PARAMR //DIV/V,N,BG0w3/C, Y,G=0.0/C,Y, W3*0,0 S 

PARAMR _//COMPLEX//RGOw3//_V f N # Gnw3 % 

PARAMR />blV/V,N,RdOwa/| ,b/c, y, waso.o s 
PARAMR //COMPLEX//ROOwa//V,N # QOW« S 
ADD B2PP , BPP/BX 1 S 
ADD 'KPP,BXT/BX2 7 g'063' S 
add KaPP,gx2/Bx3/nofia $ 

ADD BCPP,BJC3/Bxa J _ 

■ CHKPNT BX4 s 
EOUIV Bxa,B2Pp/ALwAVS S 

CHKPNT_B2PP S __ 

PARAm'l 82PP//PPESEMCE////V # N,nQB2PP S 
S generate stiffness terms 

ADD _K 2PP , KUPP/KX 1 /JLUi 0 ,.} \ «J _S 

ADD KCPP,KX1/KX2 § 

ADD Kx2, AER0/KX3//C,Y,MAEpas(i ,,o,J S 
CMKPNT K X 3 S 
EOUIV KX3,K2PP/ALWAYS S 
CHKPNT K2PP S 


PARAML K2PP//pRESENeE/ ///V ,N t NQK2PP J 

AOB'COUPLX#M2PP/M2PPX s 
CHKPNT M2PPX s 

EOUIV M 2 P P X , M 2_P P_/ A L w A Y 3 S 

CHKPNT M2PP S 

PARAML M2PP//PRESENCE////v-,N,N0M2PP S 

ALTER 106 S 

CHKPNT M2PP,B2PP,K2PP' t' 

ALTER 108,108 S 

CHKPNT M20D f 82DD,K2DD S 

ALTER 1 OR , 1 OR S 

CKAD USE TD, GM , GO, , , , , K2Pp , M2PP , B2PP/K0 D , BDD, mDO , GMO , 

G0D,K2D0 f M2DD, 8200/0^, CMPLEV/C,N*D ISP/C, N, DIRECT/C, Y»G*0. 0/C. 

N,O t O/b,0/V ( N,NOK2PP/V,N # NOM2PP/v«N#NOB2PP/V,N,MPCFl/V, 
N,SINGLE/V,N,0MIT/v,N,N0UE/V,N # nQK0GG/v#M,N0BGG/*1/-1 S 
ALTER 110,110 s 

CHKPNT KDD,BDD,mdo # gmD,GOD,K2DD,M2DD,0?OD’ S 
ALTER 113 


ORiarMAr 

"•wumfttr rA'ic W 

0F POOR QUALITY! 


add _KHH # /KHHx/(5aiS(i i 208,O l ) S 
ADD *" 8HH,/BHHx/(232,7l0S7,0,) S 

ALTER 114,115 
PARAM _//D0*N/0 S 

CEAD ■ KHHX". 0HHX . MHH , EED # CASE XX /PM I H ( 'CL AM A# QCE I GS /S » N, E IGVS S 

PARAM //DQwN/55000 S 


ALTER 118 S GENERALIZED MASS CALCULATIONS 
HPYAD MHH, PHlH, /MpHI J 

MATMUD PHIM, , , , , /PHIHCNJ,/10 S 

MPYAQ PHIHCNJ,MPHI , /GMCX/1 % 
mathoq GMCX,,, # ,/GMCXF,/2////1..2 S 
M A i 3 R N GHCXF // S 

ALTER 12b S - - ... . ..... ... ... 

PURGE GPC/ALwAyS S 

MODaCC CASEXX, ClAMA,PHIh, , ,/CLAMAX,PwIHX,CASEZZ# ./CEIG J 

CHKPNT ... CLAMAX,PhIHX-I 

EGUIV PHIHX,PHrH/ALWAYS/CLAMAX,CLAMA/ALWAY8 I 
CHKPNT PHIH f CLAMA,CASEZZ S 

EQUlV .. . CASEZZ,CASEXX/ALwAYS y 

CHKPNT CASEXX S 

ALTER 133 $ STRAIN ENERGY DATA RECOVERY 
GPFDH .... CASEXX, CPHIP^KELM, KOI CT r ECT,EQEX IN, GPECT,CLAMA»QPC r - 
BGPDT,SlL»CSTM/ONRGYR,nGPFBR/pEIG/8, Y,TINY S 
OFP unrgyr,ogpfbr//s,n # carono S 

ADD .... CPHIP,/CPHIPI/(0 1 ,«1 ... 

GPFDR CASEXX, CPHIPI,KEL m # k DICT, EC T^GEXINjGPECTiCLAMAiQPC, 
BGPOT,SlL#eSTM/OKjRGYI ,OGPFBI/rEIG/C, Y,TINY S 

OFP -CNRGYl,OGPFBI//S r N,CARONO- J 

I END OF ALTER FOR COMPLEX MODES 
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Blade Pitch Distribution 


1. Static Torsional Deflection of the Heliogyro Blade 


Whenever the Heliogvro blade root is set at a prescribed angle, the 
torsional deflection (twist distribution) of the blade involves the 
numerical solution of the nonlinear differential equation: 


- M(f ) 2 T $] * 26 ■ ■« 


where c = blade chord 

T = centrifugal tension load 

Ig = distributed blade pitch inertia 

mg » applied blade pitching moment 

6 = blade pitch angle 

fi s rotational speed 
y » radial coordinate 

The difference equation used to obtain the solution of Table 1 is: 


8 n*l * 


+I 9 > sln26 „.l 

n n+1 


■» (^P) 


where Mj - “ sin 26 tlp 

' 0 
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Table 1. Static Torsional Deflection of Heliogyro Blade for R = 6250 m. 


y/R 

(c/2) 2 

(m 2 ) 

T 

(Nt) 

X 6 

(kg-m 2 ) 

Ay 

(m) 

6 

n 

(deg) 

TENPAD 
* l-cos28 

n 

1.0 


12.337 

21.33 

0 

30.00 

.500000 

.9875 

16 

29.271 

29.85 

78.125 

30.15 

.504483 

.95 

16 

79.109 

61.63 

234.375 

30.52 

.515816 

.30 

16 

143.630 

64.09 

312.5 

31.01 

.530695 

.85 

16 

205.474 

•66.41 

312.5 

31.50 

.545986 

.80 

16 

250.827 

39.08 

312.5 

32.01 

.561857 

.75 

16 

280.432 

41.14 

312.5 

32.53 

.578415 

.70 

16 

320.648 

70.59 

312.5 

33.08 

.596796 

.65 

16 

370.500 

74.38 

312.5 

33.65 

.614138 

.60 

16 

417.065 

76.04 

312.5 

34.25 

.633543 

.55 

16 

460.246 

77.79 

312.5 

34.88 

.654078 

.50 

16 

499.979 

78.96 

312.5 

35.54 

.675809 

.45 

16 

536.188 

80.22 

312.5 

36.24 

.698808 

.40. 

16 

568.783 

81.35 

312.5 

36.96 

.723159 

.35 

16 

597.693 

82.35 

312.5 

37.73 

.748953 

.30 

16 

622.859 

83.21 

312.5 

38.54 

.776291 

.25 

16 

644.237 

83.94 

312.5 

39.39 

.805280 

o 

CM 

• 

12.5 

660.934 

61.89 

312.5 

40.52 

.844322 

. 15 

9 

672.486 

39.51 

312.5 

42.15 

.900538 

. 10 

9 

679.535 

31.76 

281.25 

43.64 

.952490 

.07 

9 

583.279 

23.86 

218.75 

44.82 

.993747 

.04 

9 

685.086 

16.79 

159.64 

45.70 

1 024306 

.02783 

9 

686.169 

54.80 

84.910 

46.18 

1.041012 

.01896 

9 


0 

48.005 

46.45 

1.050662 

.01010 

9 


0 

57.97 

46.79 

1.062309 

.0012336 

9 

686.812 

206.5 

57.97 

47.15 

1.074887 


original PAGE tt> 
OP POOH quality 












AP TEND IX L 


Solar Pressure Calculation and Other Static Deformations 


1. Solar pressure at 0.25 A.U. with inertia relief for calculation of 
the blade vertical deflection in S0L 24,1: 


P 0 (1 


U) 


where y 


M ass of Center Body ♦ Flap Hinge Brace 
Total Vehicle Mass 

1 V 


(25) (P er ^ ect reflector pressure) (reflective coefficient) 


Typical values of y and p Q are 


117.041 + 32.961 
V 4047391 ’ ' °' 370933 


P 0 ■ (16) (0.9026 x lO -3 ) (0.9136 - 0.0049) 
= 1.3123 x 10 -4 Nt/m 2 


Values for the vertical loads are shown in Table 1. 
p 3 * pAyc 

Note: No solar pressure loads should be applied to the cutouts and 

reduced loads applied to reduced chord regions. 

2. Differential Chordwise Strain resulting from manufacturing deformations 
(5 x the design limit): 

Leading Edge Strain: 5 x 10' 5 

Trailing Edge Strain: -5 x 10"^ 

3. Inplane Solar Pressure Loads due to 1° of untrimmed collective pitch at 
0.25 A.U. 


p, * (p sin 0 cos 6 )Ayc 
r l ' ,r o o o' ' 


The inplane loads are shown in Table 2. 
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Table 1. Vertical Solar Pressure Loads at p 25 A.U 


118085 


118085 

118085 


236169 


236169 


236169 


236169 


236169 

210958 


185743 

185743 


210953 


236169 

236169 


236169 


236169 


236169 

236169 


188935 


141702 

110552 

039702 


y/R 

Ay 

(in) 

c 

On) 

Tip l.O 

0 

7.6289 

.9875 

187.5 

7.6289 

.95 

187.5 

7.6289 

.90 

187.5 

7.6289 

.85 

375.0 

7.6289 

.80 

375.0 

7.6289 

.75 

375.0 

7.6289 

.70 

375.0 

7.6289 

.65 

375 0 

7.6289 

.60 

375.0 

6.81' J 

.55 

375.0 

6.0 

.50 

375.0 

6.0 

.45 

375.0 

6.8145 

.40 

375.0 

7.6289 

.35 

375.0 

7.6289 

.30 

375.0 

7.6289 

.25 

375.0 

7.6289 

.20 

375.0 

7.6289 

.15 

375.0 

7.6289 

.10 

300.0 

7.6289 

.07 

225.0 

7.6289 

.04 

175.54 

7.6289 

Tru Apex .023188 

63 04 

7.6289 
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Table 2. Inplane Solar Pressure Loads at 0.25 A.U. 


y/R 

Ay 

(m) 

c 

(m) 

0 Q 

(deg) 

?! 

(Nt) 

Tip 1.0 

0. 

7.6289 

.627 

0 

.9875 

187.5 

7.6289 

.630 

.002064 

.95 

187.5 

7.6289 

.639 

.002093 

.90 

187.5 

7.6289 

.652 

.002136 

.85 

375.0 • 

7.6289 

.665 

.004357 

o 

00 

375.0 

7.6289 

.678 

.004442 

.75 

375.0 

7.6289 

.693 

.004540 

.70 

375.0 

7.6289 

.708 

.004639 

.65 

375.0 

7.6289 

.723 

.004737 

.60 

375.0 

6.8145 

.738 

.004319 

.55 

375.0 

6.0 

.755 

.003890 

.50 

375.0 

6.0 

.773 

.003983 

.45 

375.0 

6.8145 

.790 

.004623 

.40 

375.0 

7.6289 

.808 

.005294 

.35 

375.0 

7.6289 

.829 

.005431 

.30 

375.0 

7.6289 

. 851 

.005575 

.25 

375.0 

7.6289 

.873 

.005719 

o 

CM 

• 

375.0 

7.6289 

.894 

.005857 

.15 

375.0 

7.6289 

.920 

.006027 

.10 

300.0 

7.6289 

.947 

.004963 

.07 

225.0 

7.6289 

.962 

.003781 

.04 

175.54 

7.6289 

.978 

.002999 

Truss Apex .023188 

63.04 

7.6289 

.987 

.001087 
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Table 3. 

NASTRAN Solar 

Pressure Loads at 

0. 25 A.U. 

(S0LPRS) 



FORCE 

2 

10120 

1 

1 . 

0, 

o. 

,1 10OHS 


FORCE 

2 

10500 

1 

1. 

0. 

0. 

,110085 

FORCE 

2 

11000 

1 

1. 

0. 

0. 

.1 1H085 


FORCE 

2 

11500 

1 

1 . 

n. 

0, 

,236169 


FORCE 

2 

12000 

1 

1. 

o. 

0, 

.236169 

i 

FORCF 

2 

12500 

1 

1. 

0, 

0. 

.236169 


FfJMCE 

2 

1 3000 

1 

1 . 

0. 

0. 

,236169 

j 

FORCE 

? 

1 3500 

1 

1. 

0. 

0. * 

,236169 


FORCE 

2 

1U000 

1 

1. 

0. 

0. 

.210950 


FORCE 

? 

10500 

1 

1. 

o. 

0. 

,185703 


FORCE 

2 

15000 

1 

1. 

0. 

0. 

.185703 


FORCE 

? 

15500 

1 

1. 

0. 

0. 

.210R5H 


FORCF 

2 

16000 

* \ 

1. 

0. 

0. 

.236169 


FORCE 

2 

16500 

1 

1. 

0. 

0. 

.236169 


FORCE 

2 

17000 

1 

1. 

0. 

0. 

.236169 


FORCE' 

? 

17500 

1 

t. 

0. 

0. 

.236169 


force 

2 

10000 

1 

1. 

0. 

0. 

.236169 


FORCE 

2 

10500 

1 

1. 

0. 

0, 

.236169 


FORCE 

2 

1R000 

1 

1. 

0, 

0. 

.100935 


FORCE 

2 

19300 

1 

1. 

0. 

0. 

.101702 


FORCE 

2 

19600 

1 

1. 

0. 

0, 

.11053? 


FORCE 

2 

19700 

1 

1. 

0. 

0, 

,039702 


DEFORM 

5 

1001 

-.000600 

100? 

♦,000688 



deform 

5 

1013 

-,010062 

10 10 

♦,0 l a062 



deform 

5 

1051 

-.01875 

1052 

♦.01075 




*# 

= C 1 5 3 

S # 

*(50). 

S # 

*(50) , 

3. 

*1 

s* $ 


% 

GENERATE DEFORMS C1N TF AND LE 





deform 

5 

1901 

-.01125 

1902 

♦.01125 




DEFORM 

5 

1931 

-.01125 

1932 

♦,01125 



OEFohm 

5 

1961 

-,00630a 

1962 

♦.00630U 



FORCE 

2 

10120 


1. 

, 00206 a 

0. 

,118005 


FORCE 

2 

10500 

1 

1. 

,002093 

0, 

. 1 10085 


FORCE 

2 

11000 

1 

1. 

,002136 

0. 

. 1 10005 


force 

2 

1 1500 

1 

1 . 

,O0a357 

0. 

.236169 

force 

2 

12000 

1 

1 . 

.ooaaa? 

0. 

, 23o j 69 


force 

2 

12500 

1 

1. 

, 00a5ao 

0. 

,236169 


FORCE 

2 

1 3000 

1 

1. 

,000639 

0. 

,236169 

FORCE 

2 

13500 

1 

1. 

.00U73? 

0. 

.236169 


FORCE 

2 

toooo 

1 

1. 

, 0 0 a 3 l 9 

0. 

,210950 


FORCE 

2 

1U500 

1 

1 . 

,003090 

0. 

.185703 


FORCF 

? 

15000 

\ 

1 . 

,003903 

0. 

,105703 


FQRCF 

2 

15500 

1 

1. 

,00a623 

0. 

,21095? 


FORCE 

?. 

16000 

1 

1 . 

,0«32Ra 

0. 

,236169 


FORCE 

2 

16500 

1 

1 . 

, 0 05a 3 1 

0. 

.236169 

r 1BCE 

2 

1 7000 

t 

1 . 

.005575 

0. 

,236169 


FORCE 

2 

17500 

1 

1. 

,005719 

0. 

.236169 


FORCF 

2 

10000 

1 

1 . 

,005057 

0. 

,236169 


FORCE 

2 

10500 

• 

1 . 

,006027 

0. 

,2361 e° 


FORCF 

2 

19000 

1 

1. 

,000963 

0. 

.100935 

FORCE 

2 

19.100 

1 

1. 

.003/81 

0, 

,101702 


FORCE 

? 

19600 

1 

1. 

,002999 

0. 

, 1 1 0552 


FORCE 

2 

1 9 TOO 

1 

1. 

.0010H7 

0. 

,039702 
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Aerodynamic Force Coefficients 


Inplane aerodynamic force coefficients for 30° of half-P pitch and 
at zero sun angle are calculated from the following equation: 

dp x 

K 1S = * dF" c L V 


where 



K ■ washout of a^ 

Table 1 lists the inplane aerodynamic force coefficients for a, = 30° 
at 0.7R and y * 0° at 0.25 A.U. 

Inplane and vertical aerodynamic force coefficients for 5° collective 
pitch and zero sun angle at 0.25 A.U. : 

K 1S * HT c 

S^RSSSf 

where 

dp 

r x 

de " p 0 

dp z 

de" " - 3p o sin 0 o cos2e o (e c /0 o ) 

Table 2 shows inplane and vertical aerodynamic force coefficients 6 C ■ 
at 0.7R and y ■ 0* at 0.25 A.U. 
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Table 1. Inplane Aerodynamic Force 
for a, = 30° and y - 0 


023188 


(■) 

K 

K 15 

(m/rad) 

7.6289 

.833 

0 

7.6289 

.839 

.07556 

7.6289 

.859 

.07736 

7.6289 

.884 

.07961 

7.6289 

.911 

.16409 

7.6289 

.937 

.16877 

7.6289 

.969 

.17454 

7.6289 

1.000 

.18012 

7.6289 

1.032 

.18588 

6.8145 

1.065 

.17135 

6.0 

1.105 

.15654 

6.0 

1.144 

.16206 

6.8145 

1.185 

.19066 

7.6289 

1.226 

.22083 

7.6289 

1.277 

.23001 

7.6289 

1.328 

.23920 

7.6289 

1.381 

.24875 

7.6289 

* 1.434 

.25829 

7. 628° 

1.502 

.27054 

7.6289 

1.570 

.22623 

7.6289 

1.612 

.17421 

7.6289 

1.655 

.13954 

7.6289 

1.672 

.05063 J 
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Table 2. Inplane and Vertical Aerodynamic Force Coefficients 
for 6 c = 5°, y = 0° at 0.25 A.U. 




y/R 

Ay 

(m) 

c 

(m) 

e /e 

C 0 

K 15- 

(m/rad) 

"jS 

(m/rad) 

1.0 

0 

7.6289 

.886 

0 

0 

.9875 

187.5 

7.6289 

.890 

.187714 

.043183 

.95 

187.5 

7.6289 

.903 

.187714 

.043813 

.90 

187.5 

7.6289 

.921 

.187714 

.044687 

.85 

375.0 

7,6289 

.939 

.375428 

.091120 

.80 

375.0 

7.6289 

.958 

.375428 

.092564 

.75 

375.0 

7.0289 

.979 

.375428 

.095002 

.70 

375.0 

7.6289 

1.000 

.375428 

.097040 

.65 


7.6289 

1.020 

.375428 

.099078 

.60 


6.8145 

1.042 

.335350 

.090321 

.55 

375.0 

6.0 

1.067 

.295268 

.081433 

.50 

375.0 

6.0 

1.092 

.295268 

.083341 

.45 

375.0 

6.8145 

1.116 

.335350 

.096735 

.40 

375.0 

7.6289 

1.141 

.375428 

.110722 

.35 

375.0 

7.6289 

1.171 

.375428 

.113633 

.30 

375.0 

7.6289 

1.202 

.375428 

.116642 

.25 

375.0 

7.6289 

1.233 

.375428 

.119650 

.20 

375.0 

7.6289 

1.263 

.375428 

.122561 

.15 

375.0 

7.6289 

1.300 

.375428 

.126152 

.10 

300.0 

7.6289 

1.338 

.300342 

.103871 

.07 

225.0 

7.6289 

1.359 

.225256 

.079126 

.04 

175.54 

7.6289 

1.382 

.175741 

.062777 

.023188 

63.04 

7.6289 

1.394 

.063112 

.022741 


ORIGIN* 
OF POO 
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Table 3 NASTRAN Inplane Aerodynamic Force Coefficients (AER0C0EF) 


OHIO 

AERO 

0 

1 

1 

0 



$ 

* 

s 

S 

S 

s 

S 

s 

DMJG 

AE«0 

« 0t20 

5 


10120 

1 

,07556 

OMIG 

aero 

1 0*> 0 0 

5 


tosoo 

1 

,07736 

Dm I G 

aero 

11000 

5 


11000 

1 

,07961 

OMIG 

AF WO 

1 1500 

5 


11500 

1 

,16909 

Dm i g 

AERO 

12000 

5 


12000 

1 

,16877 

OMIG 

AE»0 

12500 

5 


12500 

1 

,17959 

OMIG 

AERO 

13000 

5 


13000 

1 

,18012 

DmIG 

AF.RO 

13500 

5 


13500 

1 

, 1 8588 

OMIG 

AERO 

19000 

5 


19000 
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APPENDIX M 


Table 4. NASTRAN' InpJane and Vertical Aerodynamic Force Coefficients (AER0C0EF) 
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1 9802 

19752 

19057 



CHAP 

19803 

19751 

19803 

19753 

19881 



CRAW 

19804 

19751 

19804 

t9754 

19852 



CHAR 

19851 

19751 

19851 

19801 

19856 



CRAW 

19852 

1 9751 

19052 

1980? 

19857 



C8AR 

19856 

19751 

19856 

19803 

19851 



CBAR 

19857 

19751 

19857 

19804 

19852 



PBAR 

19751 

10 


1 ,«8 




S 

RBE2S 

TO GET MOTInNS ON 

centfrline of 

PLACE 


RBE2 

1 003 

1 0000 

1? 3 

10001 

1 0002 



PBE2 

101b 

10120 

1 ?3 

10121 

10122 



PBE2 

10S3 

10500 

123 

1 0501 

10502 



= . 

*(50), 

*(500), 

m 

■ # 

*(500) 

, *(500) 

, =s S 


*(5) 








S 

GENERATES RBF. 2 

1103 THRU 

140 3 




RRE2 

1403 

14000 

123 

14001 

14002 

14003 

1 4004 

RRE2 

1 «53 

14500 

123 

14501 

1450? 



RBE2 

1503 

15000 

123 

15001 

15002 



RRE2 

1553 

15500 

123 

15501 

15502 

15503 

15504 

RBF2 

1603 

1 6000 

123 

160C1 

16002 



*# 

•(50), 

•(500), 

*# 

•(500) 

. *(500) 

, ** % 


«(5> 








S 

GENERATES R0E2 

1653 THPl. 

1903 





fwru'UM'VirvjruMfururuxi 


RBE2 

1933 

19300 

123 

1 °30 1 

193C2 


RRE2 

1963 

1 9600 

1 ?3 

1 9601 

19602 


S 

ElASTQMFHIC RF.ARImG STIFFNESS 



CBAR 

1 9991 

1 9991 

19996 

1990() 

32 


CELAS2 

1999? 

100, 

19996 

5 

1 9997 

5 

HAT 1 

11 

120.65*9 

.3 



S 

ELEMENTS 

for center 000 y 

INERTIA 


CEL AS? 

1 99911 

190, 

19996 

1 

1 9997 

1 

CELAS2 

19995 

1 90, 

19996 

2 

19097 

2 

CELAS2 

19996 

190, 

1 9906 

3 

1 9997 

3 

S 

CONSTRAINTS to 

RFMnvE singular 

DEGREES 

OF FREEDOM 

s 

ON TRUSS 

ANO Bl AOE 




SPC1 

1 

956 

19003 

19009 

15503 

15509 

spci 

1 

956 

19751 

19752 

19753 

19759 

SPCl 

1 

956 

19801 

1900? 

19003 

19509 

SPCI 

1 

956 . 

10001 

1 0002 

10121 

10122 

SPCI 

1 

956 

1 0501 

t 0*502 

11001 

11002 

*» 

*, 

mm 

• • 

*(1000), 

*(1000) 

,*(1000) 

, * ( 1 000 ) , «s 

«<7) 







s 

GENERATES SPCS 

ON LF ANO 

TE (IF 

BLADE 


SPCI 

1 

956 

1 9301 

19302 

19601 

19602 



s heliogypo blade grids 

c 




s 

blade 

CENTER LINE GRTDS 



GRID 

20000 

200 

0. 0, 


2 

GRID 

20120 

201 

0. 0, 

0. 

2 

GRID 

20500 

205 

0. 0. 

0. 

2 

*# 

*(500) 

. *(5), 

», as J 



*(16) 






S 

GENERATES GRIDS 

21000 THRU 29QO0 


GRID 

29300 

293 

0. 0. 

8, 

2 

GRID 

29600 

296 

0. 0, 

o. 

2 

S 

BLADF 

TRAILING i 

EDGE GRIDS 



GRID 

20001 

200 

3,81905 0, 

o. 

2 

GRID 

20121 

201 

3,81005 0, 

o. 

2 

GRID 

20501 

205 

3,81005 0, 

0. 

2 

*# . 

* (500 ) 

r * (5 ) » 

as % 



= (7) 



, 



S 

GENERATES GRIDS 

21001 THRU 20501 



GRIO 

?o503 

2a5 

3, 0. 

0. 

2 

GRID 

25001 

250 

3. 0. 

0. 

2 

GRID 

25503 

255 

5. 0. 

0. 

2 

GRID 

25501 

255 

3,81005 0, 

0, 

2 

= # 

* (500 ) 

. * ( 5 ) # 

a: % 



»(6) 






S 

GENERATES GRIDS 

26001 THRU 29001 



GRID 

29301 

293 

3,81005 0, 

0. 

2 

GRID 

29601 

296 

3.81005 0, 

0. 

2 

S 

BLADE 

LEADING edge GRIDS 



GRID 

20002 

200 

-3.810050, 

0 . 

2 

GRID 

20122 

201 

-3.810050. 

0 .- 

2 

GRID 

20502 

205 

-3.810050. 

0 . 

2 

«# 

*(500) 

# * (5 ) # 

<: % 



*(7) 






% 

GENERATES GRIDS 

21002 THRU 20502 


GRID 

2o5oo 

2«5 

•3. 0, 

0, 

2 

GRID 

25002 

250 

•3, 0. 

0, 

2 

GRID 

2550 A 

255 

-3. 0, 

0 . 

2 

GRID 

25502 

255 

-3.810050, 

0, 

2 

*» 

*(500) 

. * (5 ) « 

= s S 



*(6) 






S 

generates GRIDS 

2600? THRU 29002 


GRID 

29302 

293 

-3.810050, 

8, 

2 

GRID 

29602 

296 

•3,810050. 

0. 

2 

s hfliogyhu blade truss 

Flap hinge post 

AND BUCH GRIDS 

GRID 

29700 

299 ’ 

o. 173.91 

0, 

2 

GRID 

29701 

299 

3,81005 173,91 

0 . 

2 

GRID 

29702 

299 

-3.81005173,9) 

0 . 

2 

GRID 

29751 

299 

3,81005 118,51 

-1,108 

2 

GRID 

29752 

?o9 

-3.81005118,51 

-1.108 

2 

GRID 

29753 

299 

3,81005 1 1 8 , 5 1 

1.108 

2 

GRID 

29750 

299 

-3.81005118,51 

1,108 

2 

GRIO 

29801 

206 

3,81005 63,1 1 

-2.217 

2 

GRID 

29B02 

299 

-3.8100563.11 

-2,2! 7 

2 


•fSSSSSS 




- t." ■ ■ . ? — ; 



GRID 

2RH03 

200 

3, Hi <165 

63.1 1 

2.217 

? 

GRIO 

2«0O« 

299 

•3,81 U'i563 , 1 1 

2.217 

2 

GRID 

29H51 

290 

3 , 8 1 'J u 5 

7.7i 

-3.325 

2 

GRID 

29852 

299 

-3.81<J«57.71 

-3.325 

2 

CRID 

29853 

209 

«. 

7.71 

0, 

2 

grid 

298SU 

299 

0, 

7,71 

o. 

2 

GRID 

20655 

200 

-0. 

7.71 

0. 

2 

grid 

20856 

i?o 

3.H1UU5 

7.71 

3,325 

2 

GRID 

20857 

200 . 

-3.8lau57.7l 

3.325 

2 

% 

BLADE 

ATTACHMENT AND huh 

1 GRIDS 



GRID 

29990 

290 

o. 

2.15 

o. 

2 

GRID 

29996 

200 

0. 

0. 

o, ' 

2 

GRID 

29997 

209 

o. 

0. 

0. 

2 

GRID 

20098 

2 99 

0. 

0. 

o. 

2 

S 

BUDE 

ELEMFNTS 





CHAR 

2000 

1000 

20 t 20 

20000 



CBAR 

2012 


20500 

20120 



CHAR 

2050 


21000 

20500 




*(50), 

s # 

*(500) , 

* (500 ) . 

== s 


5(15) 







$ 

GENERATES bars 

2100 Thru 

1 2850 



CfJAR 

2000 


29300 

29000 



CBAR 

2030 


29o00 

20300 



CBAR 

2060 


20700 

20600 



S 

TRAILING EDGE RODS 




cnsRno 

2001 

20121 

20001 

10 

.2220-6 


COMROD 

2013 

20501 

20121 

10 

.2220-6 


CONROD 

2051 

21001 

2050! 

10 

.2220-6 


CflNRDD 

210', 

21501 

21001 

1 0 

.2702-6 


CONROD 

2151 

22001 

21501 

10 

,3730-6 


CONROD 

2201 

22501 

22001 

10 

, U678-6 


CONROD 

2251 

2300 1 

22501 

10 

,5557-6 


CONROD 

2301 

23501 

23001 

1 0 

,6378-6 


CONROD 

2351 

2U00 1 

23501 

10 

« 7 1 fl t -6 


CONROD 

2U01 

2U501 

2U001 

10 

,7805-6 


CONROD 

2U51 

25001 

2u50 3 

10 

, 8 a 9 0 -6 


CONROD 

2501 

25503 

25 A 0 1 

10 

,0076-6 


CONROD 

2551 

26001 

25501 

10 

, 96QU-6 


CONROD 

2601 

26501 

2h00 1 

1 o 

1,007-6 


CONROD 

2h5 1 

27001 

26501 

10 

1 # oa«-6 


CONROD 

2701 

27501 

2700 1 

10 

1.080-6 


CONROD 

2751 

28001 

27501 

10 

1,113-6 


CONROD 

2801 

28501 

28001 

10 

1,136-6 


CONROD 

2851 

29001 

28501 

10 

1 ,l5a-6 


29001 

10 

1 ,16U 

-6 




CONROD 

2031 

20601 

20301 

10 

1,160-6 


CONROD 

2061 

20701 

29601 

10 

1 , t "2-6 


$ 

LEADING EDGE RODS 




CONROD 

2002 

20122 

20002 

10 

,2220-6 


CONROD 

20 1 a 

20502 

20122 

1 0 

,2220-6 


CONROD 

2052 

21002 

20502 

10 

,2220-6 


CONROD 

2102 

21502 

21002 

10 

,2702-6 


CONROD 

2152 

22002 

21502 

10 

,3730-6 





CtlNROD 

22 02 

22502 

22002 

10 

,a678-6 


CONROD 

2252 

23002 

22502 

10 

,5557-6 


CtlNROD 

2302 

23502 

2300? 

10 

,6376-6 


CONROD 

2352 

2a002 

23502 • 

10 

,7iai-6 


conrod 

2«02 

2a502 

2a0O2 

1 0 

, 78 a5-6 


Cf'NRDD 

2o5? 

25002 

2a5oa 

10 

, 3U90-6 


CtlNROD 

2502 

255oa 

25002 

10 

,9076-6 


CONROD 

2552 

26002 

25502 

10 

,960a-6 


CtlNROD 

2602 

26502 

26002 

1 0 

1,007-6 


CONROD 

2652 

27002 

26502 

10 

1,008-6 


CONROD 

2702 

27502 

27002 

10 

1,080-6 


CONROD 

2752 

28002 

27502 

to 

1,113-6 


CONROD 

2802 

26502 

2600? 

10 

1 , 1 36-6 


CONROD 

2852 

29002 

26502 

10 

1 ,150-6 


CONROD 

2902 

29302 

2900? 

10 

1,160-6 


CONROD 

2932 

29602 

2930? 

1 0 

t , 1 69-6 


CONROD 

2962 

29702 

29602 

10 , 

1,172-6 


$ 

SHEAR 

PANELS 





CSHEAR 

2009 

2009 

20001 

2000? 

20122 20121 


CSHF AR 

2019 

2009 

20121 

20122 

20502 20501 


CSHEAR 

2109 

2109 

21001 

21002 

21502 21501 


PSHf AR 

2109 

30 

1,27-6 




CSHEAR 

2159 

2009 

21501 

21502 

22002 22001 


s # 

*(50 j. 

= , 

*(300) 

, * (500 ) , 

*(500), *(500) 

, ss S 

*(a) 







s 

generate shear 

PANELS 

215V Thru 

2009 


$ 

GENERA 

TE SHEAR 

PANELS 

2159 THRO 

2009 


CSHEAR 

2«59 

2009 

2a503 

?u5oa 

25002 25001 


CSHEAR 

2509 

2009 

2500 1 

25002 

25500 25503 


CSHEAR 

2559 

2009 

25501 

25502 

26002 26001 


*C5) 

*(50) , 

s # 

*(500) 

. *(500). 

**500), - (500 ) 

•* 

M 

II 

S 

GENERa 

TFS SHEAR 

PANELS 

2609 Thru 2659 


CSHEAR 

2909 

2009 

29001 

29002 

29302 29301 


CSHEAR 

2939 

2009 

29301 

29302 

2960? 29601 


CSHEAR 

2569 

2009 

29601 

29602 

2970? 29701 


s 

TRUSS 

AND BGflH 

EI.FHFNTS 



CRAR 

2970 

2970 

29700 

29701 

29850 

2 

CRAP 

2975 

2970 

2970? 

29700 

2985a 

2 

CRAR 

2993 

2993 

29655 

29850 

a* 

2 

C B A R 

299a 

2993 

2965a 

29853 

00 

2 

C B AR 

2995 

2995 

29951 

29853 

29850 

It 

CeAR 

2996 

2995 

29652 

29855 

29850 

2 

CRAR 

2997 

2995 

29853 

29856 

29850 

2 

C9AP 

299fl 

2995 

29855 

29857 

29850 

? 

CRAR 

2999 

2999 

2R9R0 

2985a 

0] 

2 

CONROD 

?9 7\ 

29701 

29751 

10 

0,7150-6 


C0NRC10 

2972 

29702 

2975? 

1 0 

0,7150-6 


CONROD 

2973 

29701 

29753 

to 

0,7150-6 


CONROD 

297a 

29702 

2975a 

to 

0,7150-6 


CONROD 

2976 

29751 

2O801 

10 

0,7150-6 


CtlNROD 

2 977 

29752 

29802 

10 

0,7159-6 


CONROO 

2978 

2975? 

29803 

10 

0,7150*6 

SSSKfflS 

CONROD 

2979 

297S« 

?980a 

to 

0,7150-6 


COSWOD 

29? l 

29801 

29851 

10 

CQNRQO 

29? ? 

29802 

29852 

1 0 

CflNRUD 

2963 

29603 

29656 

10 

conrcjd 

298 a 

29 8 09 

29857 

1 0 

conrod 

?9Ro 

29751 

29752 

1 0 

CONROD 

2987 

29753 

2975 if 

1 0 

CONROD 

2988 

29801 

29802 

1 0 

CONHOD 

?989 

29803 

29809 

1 0 

CSHCAR 

29708 

29708 

29701 

29702 

CSHEAP 

29 709 

29708 

29 701 

29702 

cshear 

29758 

29708 

29751 

29752 

CSHF AR 

29759 

29708 

29753 

?9 759 

CSHEAR 

29808 

29708 

29801 

2980? 

Cshear 

29809 

297^8 

29803 

29809 

RRAP 

3000 

29996 

29998 

5 

RBAR 

300 1 

299QH 

900001 


s 

e l t mr nt 

PROPERTIES 


PBAR 

2970 

1 0 

9,19-9 

2.589. 

PRAR 

2993 

10 

9, 1975- 

52.589. 

♦PH2993 





♦P92993A . 9 

.9 



PBAR 

2995 

10 

7.52-5 

9,73-7 

PBAR 

2999 

1 0 

2,291 -9 

1.5926 

♦PH29R9 





♦PB2999A , 3 

,3 



PRAR 

29991 

1 1 

2.291-9 

9,567- 

♦B299Q i 





♦629991 A 3 

.3 



P S H f A R 

2009 

30 

2,59-6 


PSHE ar 

29708 

20 

2 , 59-6 


1 

LUMPED 

MASSES 



CONM? 

200 

20000 

200 

2. 

C0NM2 

201 

20120 

201 

6,976 

CONM2 

205 

20500 

205 

6,667 

CC1NM2 

210 

21000 

210 

0,382 

CQNM? 

2 15 

21500 

2 1 5 

13.129 

CONM2 

220 

22000 

220 

13.299 

CONM? 

225 

22500 

2?5 

13,950 

CONM2 

230 

2 30 0 0 

230 

1 3,5«6 

CQNM2 

235 

23500 

P35 

1.3.731 

CONM2 

290 

29000 

2a 0 

13,856 

CONM2 

2u5 

29500 

205 

12.901 

CCNM2 

250 

25000 

250 

10.935 

CONM 2 

255 

25500 

255 

12,602 

CONM? 

2b0 

26000 

260 

19,250 

C0NM2 

2S5 

2h500 

265 

19,339 

CflNM? 

2 7 0 

27000 

270 

19,399 

CnNM2 

275 

27500 

275 

19,955 

CONM? 

280 

28000 

280 

19,500 

CONM2 

285 

28500 

285 

19,535 

CONM? 

290 

29000 

290 

1 1.698 

CONM? 

293 

29300 

293 

8,79? 

CONM? 

296 

29600 

296 

6,829 

C0NM2 

297 

29700 

299 

5.076 


9,7159-6 

9,7l5a-o 

9.7159- 6 
<J,7tS«-6 

9.7159- o 
u,7i 59-6 
9, 7i5a.6 

9.7159- 6 


29752 

29751 


29759 

29753 


29802 

298 01 


2 9809 

29603 


29652 

29851 


2 9857 

29856 


12396 

12396 


123956 

123956 


2,589*5 

2.589-5 


2,589-6 

9,309-6 

♦ PB29 

♦ PB29 

9,73-7 

9.73-7 


51,0193- 

51.977-7 

♦ PB29 

♦ PB29 

P, 567-7 

1 , -b 

♦ B299 
-6299 


♦ CM20 

♦ CM20 
♦C M 20‘ 

♦ CM2 1 ( 

♦ C*21' 

♦ C M 22( 

♦ CM22> 

♦ C Mj 3 ( 

♦ C m23‘ 

♦ CM2a< 

♦ CM2«‘ 
♦CM25< 

♦ C M 2 5 ‘ 

♦ C M 26 ( 

♦ C M 2h‘ 
♦C M 27f 

♦ C m ? 7 * 

♦ CM2(»( 
♦CM2B= 
♦CM?9r 
♦CM29I- 
♦CM29# 
♦CM29: 


CQNH2 

298 

29850 

299 

32.916 



♦C*29> 

CDNM2 

299 

29997 

299 

202,091 



♦C*29‘ 

♦CH200 

.1 


21.33 



.1 


♦C M 20 1 

.1 


39.333 



.1 


♦CM205 

.1 


37,11? 



,1 


♦ C M 2 1 0 

.1 


5,558 



.1 


♦CH?tS 

.1 


71,225 



• 1 


♦CH220 

.1 


73.638 



.1 


♦C*225 

.1 


75.905 



.1 


♦C M 230 

• 1 


78.02b 



• 1 


♦CM235 

.1 


80.000 



.1 


♦C M ?90 

.1 


81 .82* 



• 1 


♦CM2U5 

.1 


60.006 



.1 


♦CM250 

.1 


38,969 



.1 


♦CM255 

.1 


63.091 



.1 


♦C M 260 

• 1 


«7,679 



.1 


♦CM26? 

.1 


80,776 

« 


• 1 


♦CH270 

• 1 


09,726 



.1 


♦CM27S 

• 1 


90.531 



• 1 


♦CM280 

.1 


91,100 



• 1 



.1 


91 .701 



.1 


♦CH290 

.1 


73.653 



• 1 


♦CM293 

.1 


55,329 



.1 


♦ CH29<> 

.1 


93,212 



• 1 


♦CM297 

.1 


53.790 



,1 


♦CM29H 

.1 


219,001 



.1 


♦CM299 

.1 


.1 



.1 


S 

HLACE 

COORDINATE SYSTEM 





C0RU2R 

2 

9000 

o. 

0. 

-R. 

o. 

0. 1, ♦CR? 

♦ CR2 

1 ,732091 1 , 

o. 




• 

S 

CHARS 

TO REMOVE 

vertical 

SINGULARITIES 

ON TRUSS 


CBAR 

29751 

29751 

29751 

29701 

29856 


2 

CBAR 

29752 

29751 

2975? 

29702 

29857 


2 

CHAR 

29753 

29751 

29753 

29701 

29851 


2 

CBAR 

29 759 

29751 

29759 

29702 

29852 


2 

CBAR 

29801 

29751 

29801 

29751 

29856 


2 

CBAR 

29802 

29751 

29002 

29752 

29057- 


2 • 

CHAR 

29803 

29751 

' ->00 3 

29753 

290S1 


2 

CHAR 

29809 

29751 

, 00 a 

29759 

29052 


2 

CBAR 

29651 

29751 

29851 

29001 

29056 


2 

CBAR 

2905? 

29751 

5? 

29802 

29057 


2 

CBAR 

29856 

29751 

• 

29803 

29051 


2 

CBAR 

29057 

29751 

29057 

29809 

29852 


2 

PHAR 

29751 

10 


1 ,«8 




S 

PHE2S 

TO GET MOTIONS DM 

CENTFRl. 

TNt hf 

RLADF 


RHE? 

2003 

20000 

123 

20001 

20002 



RBE2 

2019 

20120 

123 

20121 

20122 



RHF2 

205 1 

20500 

123 

20501 

? C 3 0 2 


JJJteiNAL page & 

II U 

* C 5 0 ) , 

*(500), 

s # 

* (500 ) , 

*( k 00) 

, ss ? 

OF POOR QUALITY 

% 

GENERATES r<0E2 

2 1 03 THRIJ 

2903 




H0F.2 

290 J 

29 0 0 0 

123 

29001 

29002 



R9F ? 

2953 

29500 

123 

29501 

2 U 502 

29503 

2<i50<t 

HBE2 

2503 

25000 

1 ? 3 

25001 

2500? 



R9F2 

2553 

25500 

123 

25501 

25502 

25503 

25509 


=xaB9 ~ - 



R5F2 

2603 

26000 

123 

2600 1 

26002 


= , 

*(51 

• (50 1 , 

*(500) 

» * , 

* (500 ) , 

*(600 ) , 

• • % 

% 

GENERATES RBF.2 

2*53 THRi.1 

i 2103 



PBE2 

2133 

29300 

123 

2H01 

2130? 


RHE2 

2«63 

29600 

123 

21601 

21602 


S 

ELASTdMfmc REARING STIFFNESS 



CBAR 

21111 

29991 

29996 

29990 

92 


CE14S2 

29992 

100 , 

29996 

5 

21997 

5 

MAT2 

1 1 

120,05*9 

.3 



1 

ELEMENTS 

FOR CENTER Rf)D Y 

INERTIA 


CFLAS2 

2999U 

1UO, 

29996 

l 

29997 

1 

CEL a S2 

29995 

190. 

29906 

2 

29997 

? 

CELAS2 

29996 

190, 

29996 

3 

29997 

3 

S 

CONSTRAINTS TtJ 

remove singular 

DEGREES 

OF FREED 

% 

ON TRUSS 

AND HI AOE 




SPCl 

1 

956 

29503 

29509 

25503 

25509 

SPC1 

1 

956 

29751 

29752 

29753 

21759 

SPCl 

1 

956 

2980 1 

21802 

29803 

21809 

spc: 

1 

956 

20001 

20002 

20121 

20122 

5PC1 

1 

956 

20501 

20502 

21001 

21002 

• » 

*, 

* > 

*(10001 , 

*dooo) 

,*(1000) 

,*(1000) 

*(7) 







S 

generates SPCS 

ON LF AND 

TE OF 

blade 


SPCl 

1 

956 

29301 

21302 

21601 

21602 


APPENDIX P 



MIKEDATA.M15B1T0 


(D 


C0R029 

loo 

1 

0 . 

7500. 

C0RP2R 

101 

1 

o. 

7906,25 

CQRD2R 

105 

1 

0 . 

7125. 

C0RD2R 

110 

1 

0 . 

6750. 

C0RD2R 

115 

1 

o. 

6375, 

C0RD2R 

120 

1 

n. 

6000, 

C0RD2R 

125 

1 

o. 

5625. 

C0R02R 

130 

1 

0 . 

5250. 

C0R02R 

135 

1 

0 . 

<1875. 

C0RD2R 

1«0 


0 . 

0500, 

CGRD2R 

t«5 

i 

a. 

0125, 

CHHD2R 

150 

i 


3750. 

CDR02R 

155 

i 

o. 

3375. 

CQRD2R 

160 

l 

o. 

3000. 

C0R02R 

1.65 

i 

o. 

2625, 

C0RD2R 

170 

l 

a. 

2250. 

COR02R 

175 

l 

a. 

1A75. 

C0R02R 

1 BO 

l 

o. 

1500, 

C0R02R 

IH5 

l 

a. 

1125. 

C0R02R 

1R0 

l 

o. 

750. 

C0RD2R 

1R3 

l 

a. 

525. 

C0R02R 

1R6 

l 

o. 

300, 

C0R02R 

199 

i 

a. 

a. 

♦CR100 

300, 

7500. 

♦ ooo . 


♦CR101 

300, 

74*06.25 

♦ «oo , 


♦CR 1 05 

300, 

7125, 

♦ 4iOO . 


♦ CR1 10 

300 , 

6750 , 

♦ 4100 . 


♦ CP 1 1 5 

300, 

6375. 

♦ 4J 0 0 , 


♦CR120 

300, 

6000. 

♦ 4100. 


♦CR125 

300 , 

5625. 

♦ 000 . 


♦CR 1 30 

300 , 

5250, 

♦ 41 0 0 . 


♦CR1 55 

300, 

4*8 75 , 

♦ 4100 . 


♦CRl«o 

300, 

41500 , 

♦ 4100. 


♦CR 1 05 

300, 

a 1 25 , 

♦ 4100. 


♦CR150 

300, 

3750. 

♦ 41 0 0 . 


♦CR 1 55 

300. 

3375. 

♦ 4100. 


♦ CP 1 60 • 

300, 

3000, 

♦ 4100. 


♦CR 1 65 

300, 

2625. 

♦ aOO . 


♦CR170 

300, 

2250. 

♦ 4100 . 


♦CR1 75 

300, 

1875. 

♦ 4100. 


♦CR 1 80 

300, 

1500. 

♦ 4100 , 


♦CRIBS 

300, 

1125, 

♦ 4109, 


♦CR 1 RO 

300, 

7S0. 

♦ 4100. 


♦ C R 1 R 5 

300, 

525. 

♦ 4100 . 


♦CK1R6 

300, 

300, 

♦ 4100 . 


♦CR 1 RR 

300, 

o. 

♦ 4100 , 


GRID 

30 

1 99 

0. 

7.71 

GRID 

31 

199 

a. 

2.15 

GRIP 

32 

199 

a. 

a. 

Dm J (j 

TFNpA') 

1 0000 

5 


OMIG 

tenpad 

10120 

5 


0 M IG 

TFNPAO 

10500 

5 


DMir, 

TENPAO 

11000 

5 


OMIC 

TFNPAH 

11500 

5 


DM1G 

TETRAD 

12000 

5 



o. 

o. 

7500. 

1000, 

♦ C R 1 

o. 

c. 

7006.25 

1 000 , 

♦ C 9 1 

a. 

0, 

7125. 

1 000 . 

♦ Cti 1 

o. 

a. 

6750, 

1 000 , 

♦ C«l 

a. 

o. 

6375. 

1 000, 

♦ CR 1 

a. 

o. 

6000, 

1 000 , 

♦ CR 1 

o. 

o. 

5625, 

1 000 , 

♦ C R 1 

0, 

o. 

5250, 

1 000, 

♦ CR 1 

o. 

o. 

0875. 

1000, 

♦ CR 1 

o. 

o. 

0500, 

1000, 

♦CR 1< 

o. 

o. 

0125, 

1 000, 

♦CR 1 i 

a. 

o. 

3750, 

moo. 

♦ CR l! 

o. 

o. 

3375. 

1000, 

♦ CRl! 

0. 

o. 

3000. 

1OO0. 

♦ CRH 

a. 

0. 

2625. 

10"0. 

♦ C»H 

a. 

o. 

2250. 

1 000 , 

♦ CR11 

a. 

o. 

1875. 

1 000 , 

♦ CRH 

0. 

a, 

1500, 

t 000, 

♦ CRH 

0. 

o. 

1125. 

moo, 

♦ CRH 

o. 

o. 

750, 

iooo. 

♦ CRH 

a. 


525, 

1 000, 

♦ CRH 

a. 

o. 

300, 

1000, 

♦C«11 

a. 

a, 

a. 

moo. 

♦ CRH 

1 

1 








OF POOR QUALTV* 


1000, 


123056 


1 000, 


123056 


1000, 


123056 


mooo 

5 

a. 


10120 

5 

a. 


1 0500 

5 

0, 


11000 

5 

a. 


11500 

5 

a. 


12000 

5 

a. 



DMIG 

TF NP AO 

12500 

5 

O^TG 

TENPAO 

1 3000 

5 

OMIG 

TE^PAO 

1 3500 

5 

D M T G 

TENiPAr) 

1 4000 

5 

DM JG 

tenpad 

1 u500 

5 

Dm I G 

TF NPAO 

15000 

5 

Dm 1 g 

TF *PAD 

15500 

5 

Dm I G 

te^pad 

16000 

5 

DMIG 

TF NPAQ 

16500 

5 

DmIG 

TENPAp 

1 7000 

5 

DMIG 

TENPAp 

17500 

5 

DMIG 

tfnpad 

1 H 0 0 0 

5 

DMIG 

TFNPAO 

18500 

5 

DMIG 

TF NPAO 

19000 

5 

DMIG 

TENPAO 

1 9300 

5 

OMIG 

TE^PAD 

19600 

5 

DMIG 

TfNPAp 

19700 

5 

DMIG 

tenpad 

1985« 

5 


12500 s 

13000 5 
1 3500 S 
iaoO0 5 
1U500 5 
15000 5 
1 5500 5 
16000 5 
16500 5 
17000 5 
17500 5 
18000 5 
18500 5 
19000 5 
19300 5 
19600 5 
t°700 5 
1985" 5 


CO c ooooooooooooooo 


o 


APPENDIX Q 


MIKEDATA.M15B2T0 





criRO?w 

200 

2 

0. 

7500, 

0. 

o. 

7500. 

1 8 86 
A V t« V « 

♦CW20' 

CnRD2R 

201 

2 

o. 

7906,25 

0, 

0. 

7 9 0 , 2 5 

1 000 , 

♦ CR20 

CORD2P 

205 

2 

0. 

7125. 

0. 

0. 

7125. 

1000, 

♦CP20 1 

CGRD2R 

210 

2 

0. 

6750. 

0. 

0, 

6750. 

1 000 , 

♦CR21( 

CORD2R 

215 

2 

0. 

6375, 

0, 

0. 

8375, 

1000. 

♦ C H 2 1 *. 

CORD2R 

220 

2 

0. 

6000. 

0. 

0. 

6t PO, 

1 000 , 

♦CR221 

CURD22 

225 

2 

0. 

5625. 

0. 

0. 

5625, 

1000, 

♦ CP221 

CQR02R 

230 

2 

0. 

5250. 

0, 

0. 

5250, 

1000, 

♦ CP23( 

CORD2R 

235 

2 

0. 

9875, 

0, 

0. 

98 75, 

1 0 0 0 , 

♦ CR23! 

COR02H 

2<i0 

2 

0. 

9500 , 

0. 

0. 

95C0, 

1 000 , 

♦CR29( 

CQR02R 

295 

2 

0. 

9125, 

0, 

0, 

9125. 

1000, 

♦ CP29«. 

COPD2R 

?S0 

2 

0. 

3750. 

0, 

0. 

3750. 

loco. 

♦CH25< 

CfIRDPR 

255 

2 

0. 

3375, 

0. 

0. 

3375, 

i oi 

♦C«25‘ 

COR 02 R 

260 

2 

0. 

3000, 

0, 

0, 

3000, 

1000, 

♦CR26( 

CGRD2R 

265 

?. 

0. 

2625. 

0. 

0, 

2825, 

1000, 

♦CR26‘ 

C0R02R 

270 

2 

0. 

2250, 

0. 

0. 

2250, 

1 000 , 

♦CW27( 

CTJRD2R 

275 

2 

0. 

1875, 

0 B 

0, 

1875. 

icon. 

♦ C R 2 7 ‘ 

CORD2R 

280 

2 

0. 

1500, 

0. 

0. 

1500, 

1 000, 

♦CR28( 

CORD2R 

285 

2 

0. 

1125. 

0, 

0. 

1125, 

1000, 

♦CR28* 

COM02R 

290 

2 

0. 

750 , 

0. 

0, 

750, 

1000, 

♦CR29< 

CORD2R 

293 

2 

C. 

525 , 

c. 

0, 

525. 

1 000, 

♦CR29I 

COR02R 

296 

2 

0. 

300, 

0. 

0. 

300 , 

1000, 

♦CR29P 

C0R02R 

299 

2 

0. 

0. 

0, 

0, 

o. 

1000, 

♦CR29* 

♦CR2CG 

300, 

7500 , 

♦ 900, 







♦CR201 

300. 

7906.25 

♦ 900 . 







♦CR205 

300, 

7125. 

♦900 . 







♦ CR210 

300. 

8750. 

♦ 900, 







♦CR215 

300, 

6 $ 75 , 

♦ a 0 0 . 






j 

♦CR220 

300. 

8000 . 

♦ 900 . 







♦CR225 

300, 

5625 , 

♦ uOO . 







♦CR230 

300, 

5250. 

♦ 900 , 







♦CR235 

300, 

9875, 

♦ 9 0 0 . 







♦CR2O0 

300, 

9500. 

♦ a 0 ft . 







♦CR2U5 

300, 

9i?5. 

♦ 900 . 







♦CR250 

300, 

3750, 

♦ 900 , 







♦CR255 

300, 

3375. 

♦ 900 . 







♦CR260 

300, 

3000. 

♦ 900 . 







♦CR265 

300 , 

2825. 

♦900 , 







♦CR270 

300 , 

2250. 

♦ 900 , 







♦CR275 

300, 

1875. 

♦ aoo . 







♦CR280 

3C0, 

1500, 

♦ 900, 







♦C^?8S 

300, 

1125. 

♦ 900, 







♦ CR2*>0 

300, 

750, 

♦ 900, 







♦CR2R3 

300, 

525, 

♦900, 







♦CR2P6 

300, 

300. 

♦ 900. 







♦CR299 

300, 

0. 

♦ 900 , 







GRID 

9ft 

299 

0. 

7,71 

1000, 


123958 



GRID 

91 

299 

0. 

2.15 

1000, 


123956 



GRID 

«2 

299 

0. 

0. 

looo. 


123958 


1 

1 

OMIC 

TP NP At) 

20000 

5 


20000 

5 

o. 



Dmig 

TENPAD 

20120 

5 


20120 

5 

0, 



QHIG 

TFNPAO 

20500 

5 


205C0 

5 

o. 



OMtG 

T P. N P * Q 

2100C 

5 


21000 

5 

0 9 



OMIG 

TE^PAq 

21500 

5 


21500 

5 

0, 



OPIG 

TP NP40 

22000 

5 


22000 

5 

0. 



OMIG 

TP.NPAD 

22500 

5 


22500 

5 

0. 





Dmig 

TENPAO 

23OO0 

5 

23000 

5 

o. 

dmjg 

TE np Ar) 

23500 

5 

23500 

5 

0. 

DMIG 

TENPAD 

2'JOOO 

5 

2«000 

5 

ft. 

DMIG 

TpNPAQ 

20500 

5 

- 29500 

5 

o. 

Dmig 

TENPAD 

25*00 

5 

25000 

5 

0, 

DMIG 

tenpad 

25500 

5 

25500 

5 

o. 

DMIG 

TE NP Ao 

26000 

5 

26000 

5 

ft. 

dmig 

TPNPAO 

26500 

5 

26500 

5 

0. 

Dmig 

TENPAD 

27000 

5 

27000 

5 

o. 

Dmig 

tenpad 

27500 

5 

27500 

5 

0. 

DMIG 

tenpad 

23000 

5 

28000 

5 

o. 

dmig 

TFKPAD 

2RS00 

5 

28500 

5 

o. 

dmig 

tenpad 

29000 

5 

2 9000 

5 

o. 

dmig 

TE^pAD 

29300 

5 

29300 

5 

0. 

DMIG 

Tt^PAO 

29600 

5 

2 9603 

5 

0. 

dmig 

tenpad 

29700 

5 

29700 

5 

ft. 

DMIG 

TENPAD 

£935<i 

5 

29*59 

5 

ft. 


© 


0? POOR 




0 

APPENDIX R 

1 MIKEDATA.HUBA0D7 




% 

CHORD I NATE TRANSFORMATIONS 





CORDPP 

9000 


o. 

0. 

0. 

0. 

o, i. 

♦CW900 

♦CR9000 

1. 

0. 

o. 






GRID 

Rooooo 

9000 

0. 

0, 

0. 

9000 



S 

MATERIAL constants 






Hi T 1 

10 

120.65+9 

.3 





MAT1 

20 


,9309+9 






HAT1 

If 

30 


,9«<l+6 






s 

S 

« 

centhifucal load 

' DATA 






RFORCE* 

1 


900000 


90ofl 


.9297183963E-2 

♦ WFOF 1 

*PFORl 

O.EC 


0 , EG 


l.EO 


t 


S 

• 

MOOES 

type bulk 

DATA 






EIGR 

2 

MfilV 

o. 

’.02 

17 

1 7 

1,-9 

♦EIGH2 

♦EIGR2 

MASS 








ASET 1 

2 

10000 

12500 

15000 

19700 




aseti 

135 

10000 

10 120 

10500 

noeo 

11500 

12000 12500 


ASET1 

135 

13000 

13500 

19000 

19500 

15000 

15500 16000 


ASETl 

135 

16500 

1 7000 

17500 

10000 

18500 

19000 19300 


ASET1 

135 

196 0 0 

19700 

19059 





ASETl 

2 

20000 

22500 

25000 

29700 




ASETl 

135 

20000 

20120 

20500 

21000 

21500 

22000 22500 


ASETl 

135 

23000 

23500 

29000 

29500 

25000 

25500 26000 


ASETl 

135 

26500 

27000 

27500 

28000 

28500 

29000 29300 


ASETl 

135 

29600 

29700 

29059 





S 

CYCLIC 

modes type data 






ASETl 

t 

10000 

10500 

1 1500 

12500 

13500 

19500 15500 


ASETl 

1 

16500 

17500 

18500 

19300 

19700 

19997 


ASETl 

2 

10000 

15500 

1996? 





EIGC 

3 

HESS 

max 



1,-8 


♦ EIGC 

♦ EIGC 







9999 


EIGC 

a 

DET 

MAX 



1.-3 


♦EIGCD 

♦E1GCD 

o. 

.02 

o. 

'.09 

.002 

3 

9 


SPC1 

l 

<156 

1 q997 






SPC1 

1 

956 

29997 






SPCl 

l 

356 

900000 






SPCl 

1 

356 

900001 






DMI 

OMI 

0 

1 

2 

? 


6 6 


OMI* 

OMI 


1 


2 


-.59000000000-1 

♦ 01 

♦ 01 









0«I* 

OMI 


2 


1 


♦,59000000000-1 

♦ 02 

♦ 02 









OMI 

0M2 

0 

6 

2 

2 


6 6 


OMI* 

0M2 


1 


1 


-.72900000000-3 

♦ Ml 

♦ hi 









Dm I « 

0M2 


2 


2 


-, 72900000000-3 

♦ M2 

DM I * 

0M2 


? 


2 


o 'do 

♦ M2 

♦ m2 









DMI* 

Om2 


5 


5 


♦.72900000000-3 

♦ M3 

♦ M 3 









OMIC 

COUPLx 

0 

1 

3 

0 







i 




~ sa ms ■ -» 1 • 


.a 

P 

. jE 

d: 

v 






LI 

1 

r.V' 

tv 

v 1 

V *** 

V * 


J *' 




■ 


M 


© 


UMIG 

cnuPLx 

900000 

1 


904490 

0 

o. 

l . 

OMIG 

CQUPLx 

9 0 1) 0 0 0 

? 


994490 

0 

1 . 

e. 

Of IG 

r.oup Lx 

949990 

0 


900000 

1 

n . 

- 1 , 


Roooon 

2 

♦ 1 . 

0. 





Dm I G 

TENPAO 

0 

6 

1 

0 




SPOIMT 

999990 








TMQ BEADED “ODEL DATA 







LU*EW 

HUB GRID 







GRID 

900001 

9000 

o. 

0 , 

-p. 

9000 




STRUCTURE BETWEEN T«F 

Twn HIJB s 





C6AR 

9 0 0 0 0 

90000 

900001 

900000 

1 000, 

o. 

0. 

l 

PBAH 

90000 

10 

3,177- 

a ,fl93-a 

,H93-« 

,893-0 




CYCLIC 

MODES DATA FCIP 

two beaded 

MODELS 



• 

Dm I G 

cnuPix 

900001 

1 


999991 

0 

o. 

l , 

dmig 

C11UPLX 

900001 

2 


999991 

0 

1, 

o. 

DMIG 

CC 1 ' °LX 

999991 

0 


900001 

1 

o. 

-l. 


♦CGuPl 


♦ COUPl 


($) 




